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tectonic  regions  of  the  wlstern  United  States.  A  variety  of  measurements  were  used  in 
the  0.3-1  Hz  intermediate  frequency  band  including  S-wave  periods,  P  and  S  rise  times, 
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Part  Tl:  <-  Analyses  of  Lon^  Period  Data  1  ~p 

Jl'ong  period  multiple  S  and  ScS  phases  observed  in  northern  Europe  were  analyzed  to 
determine  mantle  attenuation  in  the  0.02  to  0.2  Hz  range  under  the  Eurasian  shield.  Pwo 
groups  of  events  are  used:  deep  Far-Eastern  earthquakes  and  large  earthquakes  near  the 
edges  of  the  shield  areas  of  Eurasia.^  The  Q  of  the  upper  mantle  under  the  Eurasian  shield 
roeion  was  estimated  in  the  time  domain  by  taking  amplitude  ratios  and  in  the  ltquency 
domain  by  taking  spectral  amplitude  , ratios  among  various  arrivals  and  matching  waveforms 
of  synthetic  seismograms  and  observed  data.  Under  shield  regions  SS/S  amplitude  ratios 
give  t*  ro  2.5-3  seconds  and  multiple  ScS  amplitude  ratios  give  t*  ^  4.  secon  s.  nt er 
tectonic  regions,  multiple  S  amplitude  ratios  suggest  t*  'w  5  seconds.  The  results  show 
that  the  upper  mantle  Q  under  this  shield  area  is  larger  than  the  global  average,  but  less 
than  the  Q  values  inferred  from  our  studies  of  short  period  data  in  the  same  area.  Ire 
liminary  result/’ also  suggest  that  at  frequencies  around  0.02  to  0.2  Hz,  there  is  a  A  tR 
dif f erential ytfr  around  2  seconds  between  shield  and  tectonic  regions. 
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Part  III:  The  Q  Model t 

A  large  set  of  broad  band  data  was  analyzed  to  determine  the  frequency  and  depth  de¬ 
pendence  of  Q  for  P  and  S  waves  under  the  northern  shield  areas  of  Eurasia.  A  wide  range 
of  techniques  utilizing  spectra,  amplitude  ratios  and  waveform  modeling  were  used  to 
derive  apparent  and  absolute  t*  estimated  for  P  and  S  waves  covering  the  seismic  band 
between  0.01  to  10  Hz.  >  A  Q  model  of  the  Eurasian  shield  was  constructed  on  the  basis  of 
these  results.  The  dat|a  require  a  model  in  which  Q  increases  with _ frequency  and  which  is 
c^’-acteri  zed  by  Q  val/ues  in  the  upper  mantle  that  are  generally  higher  than  those  o 
global  average  models7.  The  model  with  the  best  fit  includes  a  minimum  in  Q  between  about 
100  and  200  km  depth  and  high  Q  values  of  the  order  of  thousands  throughout  the  bulk  of 
the  mantle.  The  long  period  multiple  ScS  may  require  a  a  low  0  zone  near  the  core  mant.e 
boundary.  Preliminary  results  suggest  that  t*  versus  frequency  in  tectonic  regions  is 
higher  and  roughly  parallel,  or  slightly  divergent  towards  low  frequencies  when  compared 
to  t*  versus  frequency  in  shield  regions,  with  a  t*  differential  of  about  0.2  seconds. 
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Supplement:  Methodologies  for  Estimating  t*(f)  from  Short  Period  Body  Waves  and  Regional 
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Variations  of  t*(f)  in  the  United  States. 


In  this  paper  we  discuss  some  aspects  of  estimating  t*  from  short  period  body  waves 
and  present  some  limits  on  t*(f)  models  for  the  central  and  southwestern  United  States 
(CUS  and  SWUS).  We  find  that  for  short  period  data,  with  frequencies  above  or  2 
Hz,  while  the  average  spectral  shape  is  stable,  the  smaller  details  of  the  spectra  are 
not;  thus,  only  an  average  t*,  and  not  a  frequency-dependent  t*,  can  be  derived  from  such 
information.  Also,  amplitudes  are  extremely  variable  for  short  period  data,  and  thus  a 
great  deal  of  data  from  many  stations  and  azimuths  must  be  used  when  amplitudes  are  in¬ 
cluded  in  attenuation  studies. 

The  predictions  of  three  pairs  of  models  for  t*(f)  in  the  central  and  southwestern 
United  States  are  compared  with  time  domain  observations  of  amplitudes  and  waveforms  and 
frequency  domain  observations  of  spectral  slopes  to  put  bounds  on  the  attenuation  under 
the  different  parts  of  the  country.  A  model  with  the  t*  values  of  the  CUS  an  .  con 
verging  at  low  frequencies  and  differing  slightly  at  high  frequencies  matches  the  spectral 
domain  characteristics,  but  not  the  time  domain  amplitudes  and  waveforms  of  s  ort  perioc 
body  waves.  A  model  with  t*  curves  converging  at  low  frequencies,  but  diverging  strong  y 
at  high  frequencies  matches  the  time  domain  observations,  but  not  the  spectra  :  • 

model  with  near ly-paral lei  t*(f)  curves  for  the  central  and  southwestern  United  States 
satisfies  both  the  time  and  frequency  domain  observations. 
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We  conclude  that  use  of  both  time  and  frequency  domain  information  is  essential  in 
determining  t*(f)  models.  For  the  central  and  southwestern  United  States,  a  model  with 
nearly-paral  1  el  t*(f)  curves,  where  At*  ^  0.2  seconds,  satisfies  both  kinds  oT  data 
in  the  Q. 3-2  Hz  frequency  range. 


Introduction 


In  order  to  determine  the  functional  forms  of  frequency  dependence  of  Q  in  the 
mantle  under  a  specific  region,  broad-band  measurements  of  attenuation  must  d- 
made  on  seismic  waves  that  traversed  the  mantle  under  the  same  region.  Since  the 
various  kinds  of  methods  for  the  estimation  of  Q  are  biased  in  different  ways  by  other 
factors  it  is  desirable  to  combine  various  methodologies  rather  than  to  rely  on  one 
method.  Only  Q  models  that  are  compatible  all  kinds  of  data  can  be  accepted  as  valid 
representations  of  the  Earth’s  structure.  While  the  frequency  dependence  of  Q  is  often 
asserted  on  the  incompatibility  of  Q  values  from  free  oscillation  studies  and  short 
period  body  wave  measurements,  it  must  be  remembered  that  while  free  oscillations 
give  globally  averaged  Q  structures,  attenuation  estimates  from  body  waves  are 
specific  to  given  paths.  The  discrepancies  in  the  Q  between  the  two  kinds  of  methods 
can  be  explained  in  other  ways.  This  study  attempts  to  determine  the  frequency 
dependence  under  the  Eurasian  shield  from  broad-band  analyses  of  teleseismic  body 
waves  confined  to  the  shield  area.  The  general  plan  of  the  study  is  to  make  Q  estimates 
in  several  overlapping  narrow  frequency  bands  and  piece  together  the  best  fitting 
model,  frequency  dependent  if  necessary,  from  a  variety  of  measurements. 

In  this  paper,  we  present  results  of  a  study  of  t*  under  the  Eurasian  shield  for  fre¬ 
quencies  between  0.5  and  8  Hz.  The  following  paper  presents  results  for  frequencies 
between  0.3  and  0.5  Hz  (Lees  et  al,  1984,  hereafter  referred  to  as  Paper  II)  and  the 
third  paper  presents  the  frequency  dependent  Q  models  that  we  have  developed  from 
this  broadband  study  (Der  et  al,  1984a,  hereafter  referred  to  as  Paper  111). 

Estimation  of  Q  in  the  short  period  band  is  facilitated  by  the  increasing  availability 
of  high  quality  digital  data  which  enables  one  to  measure  the  high  frequency  (3-8  Hz) 
content  of  teleseismic  short  period  body  waves.  The  sensitivity  of  short  period  body 
wave  spectra  to  Q  is  such  that  very  accurate  estimates  (to  0.1  seconds)  of  the  attenua¬ 
tion  parameter  tp  may  be  obtained  in  spite  of  uncertainties  in  source  spectra.  Even 


crude  measurements  of  some  obvious  time  domain  signal  characteristics  such  as  rise 
times  (Stewart,  19B4)  or  dominant  periods  (Der  et  al,  1982)  put  stringent  limits  on  the 
possible  values  of  mantle  Q.  The  mere  observability  of  some  short  period  arrivals,  snah 
as  SS,  is  proof  of  high  Q  along  the  path  (for  example,  SS  at  80°  with  a  period  of  less  than 
2  seconds  and  £55  ~  5  seconds  or  £5  ~  2.5  seconds  suggest  an  average  Qs  ~  350  for  the 
upper  900  km  of  the  mantle).  Clearly,  the  formerly  accepted  values  of  tp  of  1  second 
and  £5  of  4  seconds  are  incompatible  with  the  time  and  frequency  domain  characteris¬ 
tics  of  seismic  waves  in  the  intermediate  and  short  period  bands. 

The  northern  parts  of  the  Eurasian  continent  constitute  a  classical  shield  with  high 
upper  mantle  velocities.  Although  recent  work  in  tomographic  inversions  for  laterally 
varying  earth  structure  (Clayton  and  Comer,  1983;  Nataf  et  al,  1984;  Woodhouse  and 
Dziewonski,  1984)  are  rapidly  changing  our  concepts  of  tectonic  classifications,  the 
results  do  show  that  this  part  of  the  world  is  underlain  by  a  deep  high  velocity  "root". 
In  contrast,  the  adjoining  regions  of  Iran,  Afghanistan,  Korea,  China,  and  Southeast  Asia 
are  characterized  by  low  velocities  in  the  underlying  upper  mantle.  The  shield  areas  of 
this  continent  are  also  bounded  by  the  North  American  Plate  east  of  the  Lena  river. 
The  decompositions  of  the  Earth  velocity  structure  using  low  order  spherical  harmon¬ 
ics  do  not  allow  us  to  distinguish  fine  details,  and  rrliany  apparent  fine  structures  may 
be  side  lobes  produced  by  the  analysis  procedures  (Tanimoto  and  Kanamori,  1984). 
Since  we  want  to  determine  the  Q  model  for  the  mantle  under  the  shield  area,  we  shall 

take  special  care  not  to  mix  data  from  various  regions,  but  separate  out  the  effects  of 

\ 

propagation  through  the  mantle  under  the  adjoining  tectonically  active  regions. 

Previous  work  on  attenuation  in  the  short  period  band  has  yielded  tp  around  0.15 
seconds  (Bache,  1984;  Der  et  al,  1984b)  for  Soviet  test  sites  on  the  Eurasian  shield,  and 
tp  around  0.15  seconds  and  0.35  seconds  for  the  eastern  and  western  United  States, 
respectively  (Lay  and  Helmberger,  1981;  Der  et  al,  1982;  Der  et  al,  1984b). 


initial  P  wave,  a  Parzen  taper,  and  10  pomt  smoothing  on  the  power  spectra.  A  noise 
window  preceding  the  arrival  was  treated  in  a  similar  manner.  The  positioning  of  the 
signal  window  was  designed  to  avoid  an  excessive  distortion  of  the  first  arrival  by  taper¬ 
ing;  the  Parzen  taper  also  has  very  small  side  lobes,  thus  preventing  the  leakage  of  the 
signal  energy  from  low  to  high  frequencies  which  may  lead  to  the  overestimation  of 
high  frequency  energy  in  the  signals.  As  typical  of  signal  spectra  from  most  USSR 
explosions  recorded  in  Fennoscandia,  these  show  high  frequency  energy  above  the 
noise  level  to  frequencies  of  6-8  Hz.  In  order  to  avoid  problems  with  the  artificial 
apparent  high  frequency  energy  generated  by  the  variable  quantization  associated  with 
gain  ranging,  all  energy  below  1%  of  the  peak  amplitude  level  is  disregarded.  Even  with 
these  conservative  assumptions  it  appears  that  the  range  oi  valid  signal  frequencies 
extends  to  6-8  Hz  in  some  teleseismic  signals  at  NORSAR.  Due  to  the  sensitivity  of  high 
frequency  energy  to  Q,  these  spectra  provide  robust  constraints  on  any  frequency 
dependent  Q  models  of  the  upper  mantle  almost  irrespective  of  the  uncertainties  in  the 
details  of  source  spectra. 

At  NORSAR  the  assumption  of  a  cube-root  scaled  source  according  to  the  von  Seg- 
gern  and  Blandford  (1972)  model  yielded  the  apparent  tp  values  from  Soviet  Peaceful 
Nuclear  Explosions  shown  plotted  on  a  map  in  Figure  2.  Over  the  1-);B  Hz  band,  t  ~  0. 1- 
0.25  seconds.  These  low  values  indicate  a  high  Q  upper  mantle  under  the  sources  as 
well  as  under  the  NORSAR  receivers.  Although  our  choice  of  source  spectral  models 
has  some  effect  on  our  tp  estimates,  other  plausible  source  models  such  as  that  pro¬ 
posed  by  Mueller  and  Murphy  (1971)  give  essentially  the  same  result  over  this  wide  fre¬ 
quency  band.  This  is  because  both  the  von  Seggern-Blandford  and  Mueller-Murphy 
models  fall  off  as  w-2  at  high  frequencies  and  the  attenuation  is  estimated  in  the  fre¬ 
quency  range  where  the  spectra  fall  off  at  the  same  rate.  However,  the  models  of 
Helmberger  and  Hadley  (1981)  and  Lay  et  al  (1984),  which  fall  off  as  W3,  would  yield 
even  lower,  sometimes  negative,  tp  values,  which  in  itself  indicates  that  such  models  do 
not  give  an  accurate  description  of  the  seismic  source.  This  is  not  surprising  since 


Calcagnile  are  superposed  on  each  figure.  Although  the  travel  time  branches  shown  in 
these  figures  may  create  the  impression  that  they  are  largely  derived  by  subjective 
interpretation,  more  detailed  sections  constructed  from  larger  sets  of  sensors,  in  the 
manner  of  King  and  Calcagnile  (1976)  (Figures  B  to  10)  show  that,  although  the  relative 
amplitudes  of  arrivals  are  variable,  the  arrivals  themselves  are  quite  distinct  and  con¬ 
tinuous  across  the  NORSAR  array. 

Inspecting  the  principal  characteristics  of  the  various  arrivals  in  the  time  domain 
the  following  major  features  are  apparent:  the  B  branch  extends  to  about  32°  in  dis¬ 
tance  and  in  the  16-20°  range  the  arrivals  on  the  A  branch  are  about  one  half  to  one 
third  the  amplitude  of  arrivals  on  the  BC  and  CD  branches.  Both  the  KCA  velocity 
model  of  King  and  Calcagnile  (1976)  and  the  K8  velocity  model  of  Given  and  Helmberger 
(1980)  were  compared  with  these  observations.  These  two  models  is  shown  in  Figure 
11;  the  PREM  model  of  Dziewonski  and  Anderson  (19B1)  was  used  for  the  lower  mantle 
in  both  cases.  Figures  12  and  13  show  sets  of  synthetic  record  sections  computed  by 
the  WKBJ  method  (Chapman,  1978;  Dey-Sarkar  and  Chapman,  1978)  for  the  KCA  and  K8 
velocity  models,  assuming  an  elastic  earth,  a  von  Seggern  and  Blandford  source  func¬ 
tion  (1972),  and  the  short  period  NORSAR  instrument  response.  The  K8  model  contains 
a  low  velocity  zone  in  the  upper  mantle  as  seen  in  Figure  11,  which  gives  loiiv  AB/CD 
amplitude  ratios,  actually  smaller  than  observed,  and  a  shortened  B  branch  terminat¬ 
ing  at  about  23°  instead  of  32°.  Inspection  of  the  data  used  by  Given  and  Helmberger 
reveals  that  the  data  for  the  Novaya  Zemlya  explosions,  which  are  included  in  the  data 
set  used  to  derive  K8,  is  quite  different  from  the  data  from  explosions  elsewhere  in  the 
USSR  in  that  the  B  branch  is  weak  or  missing  at  distances  beyond  27°.  Thus  the  K8 
model  is  inconsistent  with  observations  from  the  Eurasian  shield  and  we  must  conclude 
that  paths  involving  Novaya  Zemlya  are  not  representative  of  our  study  area.  The  KCA 
model  shows  an  extended  B  branch,  but  the  AB/CD  amplitude  ratios  around  16-20°  are 
too  large.  However,  this  amplitude  discrepancy  can  be  explained  by  a  low  Q  layer  at 
the  depth  of  about  150  km.  Thus,  the  velocity  structure  can  be  kept  identical  to  the 


King  and.  Calcagnile  model  (1976). 

Inspecting  Figures  4  to  7  reveals  that  there  is  not  much  change  in  the  patterns  of 
arrivals  as  functions  of  the  changes  in  the  frequency  bands.  This  by  itself  is  an  indica¬ 
tion  of  a  high  Q  environment  in  the  upper  mantle.  If  the  overall  envelopes  of  all  arrivals 
were  the  same  in  all  frequency  bands  this  would  indicate  an  infinite  Q  throughout  the 
upper  mantle  under  the  Eurasian  shield. 

To  measure  more  subtle  changes  in  the  frequency  content  of  the  various  travel 
time  branches  and  to  put  further  constraints  on  the  upper  mantle  Q  model,  we  have 
computed  spectral  ratios  between  some  of  the  P  travel  time  branches.  This  approach 
has  the  advantage  that  source  and  receiver  effects  cancel  and  using  many  traces  one 
can  assess  the  reliability  of  the  results  by  comparing  the  spectral  ratios  between 
neighboring  sensors.  First,  we  need  to  consider  the  effect  of  upper  mantle  structure 
on  the  P-wave  spectra,  as  any  structural  contribution  to  the  spectra  must  be  separated 
from  the  contributions  of  intrinsic  attenuation.  To  investigate  this  problem,  avc  com¬ 
puted  suites  of  synthetic  seismograms  using  the  WKBJ  method  for  the  King  and  Calcag¬ 
nile  velocity  model.  Since  no  attenuation  was  included  in  producing  the  synthetics, 

spectral  ratios  between  branches  of  the  synthetics  give  an  estimate  of  the  structural 
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contribution  to  the  total  measured  attenuation.  Spectral  ratios  at  16° 1  and  30°  are 
shown  in  Figure  14.  For  CD/AB  at  16°,  P"  ~  -0.02  seconds  and  for  AB/EF  at  30°,  P  ~ 
0.04  seconds;  thus  at  these  high  frequencies,  upper  mantle  structure  may  have  an 

observable  effect  on  the  spectra  and  should  be  considered  when  estimating  the  intrin- 

\ 

sic  attenuation.  It  may  be  also  argued  however,  that  these  synthetic  records  may  not 
be  realistic  and  detailed  enough  to  model  fine  spectral  details.  Indeed  we  have  only 
included  the  directly  returned  rays  in  the  calculations,  and  reverberations,  which 
cause  long  ringing  wavetrains  were  not  modeled.  Since  the  synthetics  for  the  directly 
returned  rays  indicate  spectral  differences  between  the  branches  similar  to  those 
observed,  as  shown  below,  the  limiting  assumption  that  may  be  made  is  that  all  these 
differences  are  due  to  anelasticity.  This  would  overstate  the  amount  of  attenuation  in 


results  are  discussed  in  the  following  sections. 


Rise  Times  of  P  Waves  from  Nuclear  Explosions 


A  time  domain  approach  to  the  estimation  of  Q  for  body  waves  is  the  maLoh:::g  ei 
rise  times  (Stewart,  1984).  The  idea  is  that  t*  is  the  controlling  factor  that  determines 
the  overall  gross  shape  of  the  initial  swing  of  the  short  period  P  waveform  as  character¬ 
ized  by  a  "rise  time".  Other  factors  that  may  influence  the  rise  time  are  the  shape  of 
the  initial  source  pulse,  the  time  lag  between  P  and  pP,  and  "stochastic  dispersion" 
(Richards  and  Menke,  1983;  McLaughlin  et  al,  1984).  For  nuclear  explosions  we  have  a 
reasonably  good  idea  about  the  minimum  rise  times  of  the  source  pulses  from  near 
field  measurements;  these  are  the  shortest  for  shots  in  hard  rock.  The  minimum  P-pP 
time  lag  may  be  obtained  by  assuming  that  the  explosion  was  buried  at  a  relatively 
shallow  scaled  depth  (h=0.07Y1/3,  where  the  depth  h  is  given  in  krn  and  the  explosion 
yield  Y  in  kilotons)  and  assuming  that  the  near  surface  velocities  are  high.  Deepm 
burial  depths  and  lower  uphole  velocities  would  give  overly  high,  conservative  upper 
limits  for  tp. 


Rise  time  is  usually  defined  as  the  time  between  the  first  maximum  and  the  inter¬ 
section  of  the  tangent  to  the  rising  portion  of  jthe  P  waveform  with  the  maximum  slope 
with  the  zero  amplitude  level.  Since  the  maximum  slope  is  hard  to  measure  we  have 
used  the  time  between  the  first  break  and  the  first  maximum  as  rise  time.  We  have 
used  the  same  definition  in  our  theoretical  simulations  and,  therefore,  this 


modification  of  the  procedure  has  no  effect  on  out’  conclusions.  We  have  measured  the 
rise  times  for  a  set  of  USSR  nuclear  explosions  and  one  Indian  nuclear  explosion  at 
NORSAR.  The  theoretical  rise  times  were  computed  by  using  the  cube-root  scaled  von 
Seggern  andBlandford  (1972)  granite  source  model,  the  minimum  scale  depth,  a  range 
of  apparent  t *  from  0  to  1  seconds,  a  surface  reflection  coefficient  of  unity  and  an 
uphole  velocity  of  5.5  km/sec.  The  observed  rise  times  are  plotted  against  the 
apparent  t*  values  in  Figure  18.  The  four  curves  in  the  figure  correspond  to  the 


attributing  the  entire  difference  in  period  between  the  SS  and  S  phases  to  attenuation, 
the  estimate  of  P  is  an  upper  bound. 

As  above,  comparison  of  the  measured  periods  with  the  periods  of  synthetic  simu 
lations  such  as  those  in  Figure  20  and  Table  1  were  used  to  estimate  £5  for  each  of  the  S 
and  SS.  Consideration  of  different  source  durations  for  the  synthetic  waveforms  puts 
bounds  on  the  estimates  of  A £5  for  each  S-SS  pair.  The  changes  in  period  of  around  ton 
S-SS  pairs  that  we  have  observed  give  differentials  of  less  them  3  seconds  for  S  and 
SS  with  frequencies  in  the  0.3-1  Hz  range.  This  differential  measurement  gives  the 
additional  attenuation  of  the  SS  relative  to  the  S  from  the  passage  of  SS  through  the 
upper  mantle  to  its  surface  reflection  point  and  back  down  through  the  upper  mantle. 
Thus  the  measurement  corresponds  to  the  ts  for  a  double  pass  through  the  mantle  in 
the  region  of  the  SS  surface  reflection  point. 

Rise  Times  of  S  Waves  from  Earthquakes 

The  rise  times  of  the  short  period  S  arrivals  from  deep  events  can  also  be  used  to 
give  an  estimate  of  £5.  As  discussed  earlier  with  relation  to  short  period  P  waves,  t* 
measurements  from  rise  times  give  upper  bounds  on  estimates  of  t*  due  to  scattering 
of  high  frequency  energy  from  the  initial  pulse  into  the  fcoda  and  the  resultant 
broadening  of  the  initial  part  of  the  waveform. 

We  again  estimate  £j  by  comparison  of  the  observed  pulses  with  synthetic 
waveforms,  such  as  tnose  shown  in  Figure  20,  generated  for  a  variety  of  source  dura- 
tions  and  £  s.  Table  11  is  a  tabulation  of  the  rise  times  of  such  synthetics.  Figure  22  is 
a  plot  of  £*  versus  rise  time  for  a  range  of  source  durations.  The  arrows  on  the  right 
hand  side  of  the  figure  correspond  to  the  rise  times  measured  from  deep  Far  Eastern 
events  as  observed  at  Northern  European  and  Scandinavian  stations.  For  source  dura- 
tions  of  1  or  2  seconds,  these  short  periods  cannot  be  reproduced  for  £5  much  larger 
than  1  to  2  seconds.  Doubling  this  value  one  gets  an  upper  limit  for  £5  of  2-4  seconds 


20  and  Table  I,  and  indicated  that  the  associated  tj>  differential  must  be  less  than  0.5 
seconds.  This  is  in  good  agreement  with  the  findings  given  in  the  previous  section. 


Nine  sets  of  measurements  gave  A ts  <  2.5  seconds  for  a  double  pass  through  the  upper 
mantle. 

Comparison  of  SeS-ScP  Pairs 

Following  the  approach  by  Burdick  (1983),  we  also  looked  for  ScS-ScP  pairs  on 
short  period  records  of  deep  earthquakes.  Several  pairs  of  the  observed  ScS-ScP  pairs 
are  shown  in  Figure  25.  While  the  paths  of  the  two  phases  have  different  reflection 
points  on  the  core-mantle  boundary,  they  follow  nearly  the  same  path  in  the  upper 
mantle  source  region  where  more  lateral  heterogeneities  are  expected  to  exist. 
Because  the  first  leg  of  both  phases  is  an  S  wave,  the  difference  in  t*  between  ScS  and 
ScP  is  simply  due  to  the  difference  in  attenuation  between  the  second  legs  of  the  two 
phases.  Thus,  this  measurement  is  mostly  sensitive  to  attenuation  in  the  upper1  mantle 
in  the  vicinity  of  the  receiver,  and  is  relatively  independent  of  the  nature  of  the  source 
region. 

The  differences  in  period  of  the  ScS  and  ScP  observations  were  used  to  estimate 
tj.  For  each  phase,  t T  was  estimated  from  the  periods  of  synthetics,  like  those  in  Fig¬ 
ure  20  and  Table  I,  generated  for  a  range  of  sourdte  durations  and  t* s.  Assuming  all 
losses  in  compression,  t*$  -  4 tp.  If  ts  is  the  t*  for  an  S  wave  traveling  from  the  core 
mantle  boundary  to  the  surface,  £s(ScS)  =  2£s  and  £s(ScP)  =  t$  +  tp  =  (5/4 )t$.  There¬ 
fore,  A V  =  F(ScS)  -  P(ScP)  =  2tJ  -  (5/4 )fj  =  (3/4 )£J.  For  a  single  pass  through  the 
upper  mantle  (core  to  surface),  ts  =  (4/3)A£ ',  and  for  a  double  pass  through  the  upper 
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mantle  (say  ScS  from  a  surface  source),  ts  =  (8/3)A£  . 

Measurements  from  several  ScS-ScP  pairs  give  ts  ~  2.5  seconds  for  a  double  pass 
through  the  mantle  at  around  0.5  Hz.  Our  measurements  are  generally  smaller  than 
Burdick’s  result  of  tj  ~  4  to  7  seconds  because  the  differences  in  the  dominant  periods 
of  ScS  and  ScP  phases  are  less  than  in  Burdick's  study.  Since  our  observed  ScS  and 


FcP  phases  tended  to  be  very  small,  this  measurement  alone  is  not  sufficient  to  define 
t*  in  the  mid-period  band,  though  these  results  are  consistent  with  our  other  ones  for 
the  same  frequency  range. 

Summary  and  Conclusions 

A  variety  of  different  techniques  were  used  to  estimate  t*  under  the  Eurasian 
shield  in  the  frequency  range  from  0.3  to  8  Hz.  At  the  highest  frequencies,  recordings 
from  NORSAR  of  Soviet  nuclear  explosions  were  used.  Spectra  give  tp  ~  0.15  seconds. 
Amplitudes  and  spectra  of  the  branches  of  suites  of  record  sections  suggest  that  a  low- 
velocity  zone  is  not  required  in  the  upper  mantle  beneath  the  Eurasian  shield,  while  a 
slight  low  Q  layer  is  needed  to  explain  the  branch  amplitude  ratios.  The  rise  times  of  I3 
waves  from -m-olear  explosions  give  higher  estimates  of  tp,  around  0.5  seconds,  but  this 
may  be  due  to  scattering  of  the  initial  high  energy  in  the  P-wave  train  back  into  the 
coda.  The  periods  and  rise  times  of  S  waves  and  differential  periods  of  P-PP,  S-SS,  and 
ScP-ScS  pairs  all  give  f/  ~  2  to  3  seconds  for  frequencies  around  0.3  to  1  Hz.  And,  com¬ 
parison  of  short  and  long  period  SS/S  amplitude  ratios  gives  tj  around  2  to  2.5  seconds 
across  the  frequency  range  of  0.07  to  0.5  Hz.  |i 

Using  the  assumption  t$  =  Atp  to  relate  the  tp  and  t$  measurements,  it  is  clear 
that  there  is  a  definite  difference  in  t*  above  1  Hz  and  below  1  Hz,  although  the  data 
cannot  resolve  the  details  of  the  transition  between  t*  values.  Furthermore,  all  of  the 
t*  values  estimated  from  the  data  in  this  study  are  significantly  below  the  whole  earth 
estimates  at  long  periods  from  free  oscillations  (tp  =  1  second  and  t$  =  4  seconds),  so 
the  frequency  dependence  in  t  *  extends  to  lower  frequencies. 
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Figure  Captions 


Figure  1.  Representative  P-wave  spectra  at  NORSAR  from  Soviet  peaceful  nuclear' 
explosions. 


Figure  2.  Locations  of  Soviet  peaceful  nuclear  explosions  used  in  this  study.  Each  PNE 
is  labeled  with  the  value  of  tp  for  the  path  to  NORSAR,  assuming  a  von  Seggern  and 
Blandford  source  model. 


Figure  3.  Array  averaged  power  spectra  at  NORSAR  of  (a)  a  Soviet  peaceful  nuclear 
explosion  on  24  November  1972  and  (b)  a  Kazakh  explosion  on  10  July  19  ^3.  The 
bounds  are  95%  confidence  limits.  Theoretical  spectra  derived  using  a  von  Seggern  and 
Blandford  (1972)  source  model  and  the  t  *  noted  nest  to  each  figure. 


Figure  4.  Seismic  profiles  from  Soviet  PNE's  at  NORSAR.  The  travel  time  triplications 
from  the  model  of  King  and  Calcagnile  (1976)  are  superposed. 


Figure  5.  Seismic  profiles  from  Soviet  PNE’s  at  NORSAR,  filtered  into  the  0.5-1. 5  Hz  fre¬ 
quency  band.  The  travel  time  triplications  from  the  model  of  King  and  Calcagnile 
(1976)  are  superposed. 

Figure  6.  Seismic  profiles  from  Soviet  PNE's  at  NORSAR,  filtered  into  the  2.0-4. 0  Hz  fre¬ 
quency  band.  The  travel  time  triplications  from  the  model  of  King  and  Calcagnile 
(1976)  are  superposed. 


Figure  7.  Seismic  profiles  from  Soviet  PNE's  at  NORSAR,  filtered  into  the  5. 0-8.0  Hz  fre¬ 
quency  band.  The  travel  time  triplications  from  the  model  of 'i  King  and  Calcagnile 
(1976)  are  superposed. 


Figure  B.  Detailed  section  of  profiles  from  Soviet  PNE's  recorded  at  NORSAR  at  around 
16°. 


Figure  9.  Detailed  section  of  profiles  from  Soviet  PNE’s  recorded  at  NORSAR  at  around 
28°. 


Figure  10.  Detailed  section  of  profiles  from  Soviet  PNE’s  recorded  at  NORSAR  at 
around  30°. 


Figure  11.  Velocity  versus  depth  for  the  KCA  (King  and  Calcagnile,  1976)  and  K8  (Given 
and  Helmberger,  1980)  P-wave  velocity  models  for  Northwestern  Eurasia. 


Figure  12.  Suite  of  synthetic  P-wave  seismograms  for  the  KCA  velocity  mode!  (King  end 
Calcagnile,  1976),  with  no  attenuation.  A  von  Seggern  and  Blandford  (1972)  explosion 
source  has  been  convolved  with  the  impulse  response. 

Figure  13.  Suite  of  synthetic  P-wave  seismograms  for  the  K8  velocity  model  (Given  ami 
Helmberger,  1980),  with  no  attenuation.  A  von  Seggern  and  Blandford  (1972)  explosion 
source  has  been  convolved  with  the  impulse  response. 

Figure  14.  Spectral  ratios  of  the  CD/AB  branches  at  16°  and  the  AB/EF  branches  at  o0 
for  synthetic  record  sections  generated  using  the  KCA  model  (King  and  Calcagnile, 
1976)  and  an  elastic  earth. 

Figure  15.  Representative  spectral  ratios  of  (a)  the  CD/AB  branches  at  16  and  (b)  the 
AB/EF  branches’ at  30°  for  Soviet  PNE’s  recorded  at  NORSAR.  These  spectral  ratios  give 
estimates  of  tj  of  a)  -0.07  and  b)  0.07  seconds. 

Figure  16.  Model  of  Q  versus  depth  that  satisfies  the  spectral  characteristics  of  the  P 
wave  profiles  (dashed  line)  as  compared  to  a  model  of  the  western  United  States  (Dei 
and  McElfresh,  1977). 

Figure  17.  P-wave  travel  time  curve  for  the  KCA  model  (King  and  Calcagnile,  1970)  >v“hh 
tp  values  labeled  for  various  points  on  the  travel  time  curve.  The  Q  model  in  Figure  15 
was  used  in  a  raytracing  program  to  calculate  the  tp  values. 

Figure  IB.  Estimation  of  tj>  from  the  rise  times  of  Soviet  PNE's  as  observed  at  NORSAR. 
The  lines  are  theoretical  curves  of  rise  time  versus  tj  for  different  yields,  derived  from 

Measured  rise  times  of  synthetic  waveforms,  assuming  a  von  Seggern  and  Blandfoid 
972)  granite  source  model,  the  minimum  scale  depth,  a  range  of  t  ,  a  surface 
reflection  coefficient  of  unity,  and  an  uphole  velocity  of  5.5  km/sec. 

I 

Figure  19.  Representative  short  period  S  and  ScS  waveforms  recorded  across  Fennos- 
candia  from  events  in  the  Far  East  and  the  Hindu  Kush. 

\ 

Figure  20.  Synthetic  waveforms  for  a  range  of  source  models  and  values  of  t  *.  A  short 
period  WWSSN  instrument  response  is  convolved  with  the  waveforms  and  the  ampli¬ 
tudes  are  normalized. 

Figu  -e  21.  Estimation  of  fj  from  the  dominant  periods  of  short  period  S  waves  like 
those  in  Figure  19.  The  periods  of  synthetic  waveforms  are  plotted  on  the  graph  for  a 
range  of  source  durations  and  t  *  values.  The  observed  periods  are  marked  with  arrows 
on  the  right-hand  side  of  the  graph,  and  correspond  to  ~  1-2  seconds  for  source 
durations  between  1  and  3  seconds. 


Figure  22.  Estimation  of  ts  from  rise  Limes  of  short,  period  S  waves  like 
19.  The  rise  limes  of  syntbetie  waveforms  are  plotted  on  the  graph 
source  durations  and  t  values.  The  observed  rise  times  are  marked 
the  right-hand  side  of  the  graph,  and  correspond  to  t-s  ~  1-3  seconds. 


Figure  23.  Tracings  of  short  period  P  and  PP  arrivals  from  Far  Eastern  e 


Figure  24.  Tracings  of  short  and  long  period  S  and  SS  arrivals  from  Far 


Figure  25.  Tracings  of  ScP  and  ScS  arrivals  from  Far  Eastern  events. 
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Figure  22 


Abslracl 


Long  period  multiple  S  and  ScS  phases  observed  in  northern  Europe  were  analyzed 
to  determine  mantle  attenuation  in  the  0.02  to  0.2  Hz  range  under  the  Eurasian  shield. 
Two  groups  of  events  are  used:  deep  Far-Eastern  earthquakes  and  large  earthquakes 
near  the  edges  of  the  shield  areas  of  Eurasia.  The  Q  of  the  upper  mantle  under  the 
Eurasian  shield  region  was  estimated  in  the  time  domain  by  taking  amplitude  ratios 
and  in  the  frequency  domain  by  taking  spectral  amplitude  ratios  among  various 
arrivals  and  matching  waveforms  of  synthetic  seismograms  and  observed  data  Under 
shield  regions,  SS/S  amplitude  ratios  give  t$  ~  2.5-3  seconds  and  multiple  ScS  ampli¬ 
tude  ratios  give  4.2  seconds.  Under  tectonic  regions,  multiple  S  amplitude  ratios 
suggest  tg  ~  5  seconds.  The  results  show  that  the  upper  mantle  Q  under  this  shield 
area  is  larger  than  the  global  average,  but  less  than  the  Q  values  inferred  from  our  stu¬ 
dies  of  short  period  data  in  the  same  area.  Preliminary  results  also  suggest  that  at  fre¬ 
quencies  around  0.02  to  0.2  Hz,  there  is  a  A £5  differential  of  around  2  seconds  between 
shield  and  tectonic  regions. 


many  events.  We  have  verified  that  the  fault  plane  orientations  of  the  events  analyzed 
were  varied  enough  to  eliminate  any  significant  bias  due  to  similar  radiation  patterns, 
even  though  a  part  of  the  assumed  source  mechanisms  may  be  in  error.  We  find  that 
our  approach  is  preferable  to  only  doing  detailed  modeling  of  a  few  selected  events 
with  the  inevitable  tradeoffs  between  mantle  Q  and  source  directivity.  In  order  to 
minimize  structural  effects  we  have  considered  a  variety  of  plausible  upper  mantle 
velocity  structures.  Ideally  both  the  sources  and  receivers  should  be  located  in  shield 
areas,  but  since  practically  all  natural  earthquakes  occur  in  tectonic  areas  we  had  to 
pay  careful  attention  to  the  locations  of  surface  reflection  points  for  all  multiple  S 
phases.  We  have  also  found  that  because  of  the  complex  nature  of  long  period  S  wave 
arrivals  due  to  upper  mantle  triplications  and  scattered  background  noise,  the  spectral 
and  waveform  methods  for  estimating  Q  are  much  less  stable  than  in  the  short  period 
band.  Multiple  P  arrivals  in  the  long  period  band,  although  prominent  on  the  seismo¬ 
grams,  are  not  useful  because  of  their  relative  insensitivity  to  Qp  and  the  P-SV  conver¬ 
sions  at  reflections  which  add  the  complication  of  estimating  conversion  coefficients. 
To  simplify  our  analyses  and  reduce  structural  effects  we  have  also  restricted  our  stu¬ 
dies  to  the  SH  component  of  all  shear  wave  arrivals.  An  increase  in  £5  of  4  seconds 
should  result  in  a  change  of  S  wave  amplitudes  by  roughly  a  factor  of  2  at  the  dominant| 
periods  of  long  period  S  waves  on  the  WWSSN  records,  around  20-25  seconds. 

In  the  initial  part  of  our  work  we  have  utilized  northern  European  recordings  of 
Far-Eastern  deep  earthquakes  which  contained  many  well-separated  multiple  S  and 
ScS  arrivals.  The  amplitudes  of  the  multiple  S  arrivals  were  compared  with  full  wave 
theory  synthetic  seismograms  to  estimate  t*.  After  having  gained  some  experience 
with  the  problems  associated  with  estimating  mantle  Q  in  the  long  period  band,  we 
decided  to  analyze  a  larger  data  set  of  central  Asian  earthquakes,  mostly  not  deep,  in 
order  to  increase  the  statistical  stability  of  our  results  and  decrease  the  number  of 
tectonic-shield  mixed  paths.  The  use  of  the  relatively  uncomplicated  SS/S  amplitude 
ratios  for  these  events  tends  to  minimize  the  effects  of  structural  variations. 
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We  find  tg  ~  2.5  to  3  seconds  for  S  waves  and  4.2  seconds  for  ScS  under  shield 
regions.  Preliminary  results  suggest  that  tg  under  tectonic  regions  is  ~  5  seconds  for 
S.  These  values  of  t*  are  substantially  higher  than  those  we  found  in  the  previous 
paper  for  high  frequency  (f>l  Hz)  waves  under  the  same  region  (Dcr  ct  al,  1505a. 
hereafter  referred  to  as  Paper  1).  In  the  third  paper  of  this  series  (Der  et  al,  1985b, 
hereafter  referred  to  as  Paper  II),  we  develop  a  Q(f)  model  which  is  consistent  with 
both  the  short  and  long  period  results  for  t  * . 

Analyses  of  Data  from  Far-Eastern  Deep  Earthquakes  " 

Multiple  S  and  ScS  arrivals  from  several  deep  earthquakes  in  the  Far  East  were 
obtained  from  digital  recordings  at  northern  European  and  Scandinavian  stations.  To 
simplify  the  analysis,  the  records  were  rotated,  and  only  the  SH  components  were 
used,  which  for  the  case  of  radial  symmetry  do  not  have  conversions  to  P  waves  or  cou¬ 
pling  to  PL,  resulting  in  simpler  seismograms  and  total  reflection  at  the  core  mantle 
boundary.  The  events  are  listed  in  Table  I,  and  the  source  and  receiver  locations  are 
plotted  in  Figure  1.  Figure  1  also  shows  the  surface  reflection  points  under  the  Eura¬ 
sian  continent  of  SS  and  SSS.  The  SS  surface  reflection  points  are  clearly  under  the 
Eurasian  shield,  while  the  higher  order  multiple  S  bounce  points  are  often  under  tec¬ 
tonically  active  regions  of  China. 

Long  period  SH  seismograms  have  been  synthesized  using  full  wave  theory  for 
comparison  with  selected  rotated  records.  In  full  wave  theory,  the  synthetics  are  pro¬ 
duced  by  numerical  integration  along  the  raypath  in  the  complex  plane;  this  allows  for 
the  proper  treatment  of  such  non-ray  theory  effects  such  as  tunneling,  head  waves, 
diffracted  phases,  and  caustics  (Richards,  1973;  Choy,  1977;  Cormier  and  Richards, 
1977).  We  have  modeled  some  of  the  main  S  phases  of  the  best  four  records  from  three 
Far  Eastern  events  including  all  major  arrivals  due  to  triplications  from  reflections  off 
of  the  two  major  upper  mantle  discontinuities. 


Figure  2  shows  four  of  the  observed  seismograms  with  the  corresponding  synthetic 
seismograms.  The  synthetic  seismograms  were  computed  assuming  a  perfectly  elastic 
earth  and  using  the  SNA  shear  velocity  upper  mantle  model  of  Grand  and  Helmberger 
(1984)  merged  with  PREM  model  (Dziewonski  and  Anderson,  1981)  in  the  lower  mantle. 
The  seismograms  include  the  best  fitting  double-couple  radiation  pattern  for  the  events 
but  assume  no  directional  variations  in  the  shapes  of  the  source  time  functions,  which, 
for  these  synthetic  records  were  assumed  to  be  delta  functions.  Note  that  the  arrival 
times  of  phases  on  the  synthetics  are  close  to  those  seen  on  the  records  for  S,  ScS,  and 
SS  but  the  later  phases  are  delayed  and  their  periods  increase  on  the  observed  records 
relative  to  the  synthetics.  This  can  be  explained  by  the  observation  that  the  multiple  S 
and  ScS  arrivals  of  higher  order  cross  the  upper  mantle  in  tectonic  areas  with  lower 
velocities  (Grand  and  Helmberger,  1984)  and  lower  Q,  while  the  synthetic  records  were 
derived  from  a  high  velocity  mantle  model  for  shields.  The  synthetics  also  show  the 
complexity  of  the  SSS  waveform  in  this  distance  range  due  to  multiple  travel  paths  in 
the  upper  mantle. 

The  differences  in  peak  to  peak  amplitudes  and  pulse  shapes  of  the  synthetic 
pulses  as  compared  with  the  corresponding  observed  phases  should  thus  provide 
opportunities  to  estimate  |*.  Before  the  observed  and  synthetic  phases  can  be  com¬ 
pared,  an  appropriate  source  time  function  must  be  convolved  with  the  synthetic 
pulses.  The  source  time  function  was  modeled  by  constructing  a  transfer  function 
between  the  observed  and  synthetic  S  or  S+ScS,  depending  on  the  distance  range,  and 
convolving  the  other  synthetic' waveforms  with  the  appropriate  transfer  function.  In 
Figure  3,  some  examples  of  the  observed  SS  and  SSS  waveforms  are  shown,  matched 
with  corresponding  synthetic  seismograms  which  include  the  source  transfer  function 
and  have  been  convolved  with  attenuation  operators  corresponding  to  various  values  of 
t*.  Comparing  the  amplitudes  of  the  synthetic  and  observed  phases  gives  estimates  of 
A  t*s  4  seconds  from  SS,  but  much  higher  values,  2A£s  ~  10  seconds  or  A t$  ~  5 
seconds  from  SSS.  Thus,  it  is  clear  that  the  SSS  waves  are  being  attenuated  much 
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more  than  the  SS  waves,  as  was  noted  when  the  synthetic  records  were  compared  with 
the  observed  seismograms.  One  reason  for  this  is  that  the  SSS  have  their  first  surface 
'reflection  point  under  tectonic  regions  in  Southeast  Asia,  while  the  SS  reflection  point 
is  under  a  shield  region  as  seen  in  Figure  1.  Also,  at  this  distance  range  the  SSS  is  a 
more  complex  waveform  due  to  multiple  travel  paths  in  the  upper  mantle.  Care  must 
be  taken  to  avoid  cases  where  there  is  interference  from  sSS. 

To  complement  the  results  from  amplitude  studies,  spectral  ratios  between  multi¬ 
ple  S  or  multiple  ScS  pairs  can  be  used  to  estimate  t '  for  the  same  set  of  arrivals.  In 
order  to  more  accurately  estimate  the  intrinsic  attenuation,  any  effect  of  the  earth 
structure  on  the  spectra  must  also  be  taken  into  account.  This  is  especially  important 
with  the  SSS  and  higher  order  multiple  S  phases  since  they  interact  strongly  with  the 
upper  mantle  discontinuities  around  400  and  670  km  and  with  the  low  velocity  zone. 

The  structural  effect  on  the  spectra  was  estimated  by  taking  the  SSS/SS  spectral 
ratios  of  synthetic  pulses  for  an  elastic  earth.  In  this  case,  any  curvature  in  the  spec¬ 
tral  ratio  would  be  due  to  structural  effects.  The  synthetic  phases  were  treated  the 
same  way  as  the  data;  time  domain  phases  were  tapered  with  a  20%  cosine  window  and 
Fourier  transformed.  Three  SSS/SS  spectral  amplitude  ratios  from  synthetic  seismo- 
i grams  generated  with  £*=0  seconds  are  shown  in  Figure  4.  Within  the  range  of 
significant  instrument  response  (about  0  to  1.5  Hz),  the  resultant  spectral  ratios  are 
quite  flat,  suggesting  that  the  mantle  structure  makes  no  significant  contribution  to  t$ 
in  the  long  period  band. 

Spectral  ratios  were  taken  of  both  the  long  period  multiple  S  and  ScS  phases  from 
several  events.  The  range  of  good  signal-to-noise  ratio  varied,  but  generally  included 
0.3-0.6  Hz.  The  multiple  S  ratios  gave  a  range  of  t* s,  with  £5  averaging  around  4 
seconds  for  rays  passing  through  tectonic-region  upper  mantle.  The  ScS  spectral 
ratios  were  very  variable,  averaging  around  -2  seconds.  Representative  spectral  ratios 
are  shown  in  Figure  5.  Waveform  matching  studies  for  multiple  S  and  ScS  phases  have 


also  been  used  to  estimate  t^\  examples  of  these  fits  are  shown  in  Figure  6.  Because  of 
the  many  factors  involved  in  correcting  the  amplitudes,  such  as  geometrical  spreading, 
instrument  response,  and  radiation  pattern,  only  the  periods  were  compared  in  the 
waveform  matching.  At  least  for  short  periods,  frequency  content  tends  to  be  a  more 
robust  measurement  than  amplitude  anyway  (e.g.,  Der  and  Lees,  1984).  The  waveform 
matching  suggests  A t  *  values  for  both  S  and  ScS  that  are  low,  £  ~  1-3  seconds,  espe¬ 
cially  given  that  the  ratios  mostly  apply  to  tectonic  paths. 

The  somewhat  conflicting  results  between  the  various  approaches  used  for  Q  esti¬ 
mation  above  need  some  discussion.  The  amplitude  studies  of  SS  which  reflect  under 
the  shield  give  £5  ~  1.5  seconds,  while  SSS  which  have  a  reflection  point  under  a  tec¬ 
tonic  region  give  2 £5  ~  10  seconds  or  £5  ~  5  seconds.  These  results  are  not  surprising 
as  one  might  expect  a  difference  in  £  *  between  paths  through  shield  and  tectonic 
regions.  However,  while  the  spectral  estimates  of  multiple  S  phases  for  these  same 
paths  give  t*  ~  4  seconds  for  tectonic  paths,  the  spectral  estimates  for  the  ScS  phases 
give  t*  ~  -2  seconds.  Also,  for  the  S  phases,  spectral  ratios  of  synthetic  waveforms  for 
the  anelastic  case  (£*  =  0  seconds)  with  observed  phases  give  £  ~  0.5  seconds.  Clearly, 
some  of  the  very  low  £5  estimates  from  spectral  ratio  measurements  and  waveform 
matches  are  suspect.  It  is  likely  that  these  phases  are  either  so  complex  or  too  con¬ 
taminated  by  the  scattered  background  that  they  cannot  be  isolated  well  enough  by 
windowing  to  yield  stable,  meaningful  spectral  ratios.  This  is  suggested  by  the  fairly 
narrow  frequency  ranges  over  which  there  was  good  signal-to-noise  ratio  for  most  of 
the  arrivals.  Also,  the  long  period  instruments  are  narrowly  peaked  around  0.05  Hz,  so 
small  changes  in  period  are  not  very  discernible.  The  ScS  phases  seem  especially 
effected,  as  both  the  periods  and  amplitudes  show  a  great  deal  of  scatter;  Sipkin  and 
Jordan  (1979,1980)  have  also  observed  that  the  quality  of  the  multiple  ScS  arrivals  is 
generally  not  very  good,  perhaps  due  to  surface  scattering  or  contamination  by  surface 


wave  modes. 


Given  these  limitations  on  the  data,  it  appeared  at  this  point  of  the  investigation 
that  we  needed  a  more  robust  approach  to  define  Q  in  the  long  period  band  under 
shields. 

Analyses  of  WWSSN  Seismograms  for  Central  Asian  Earthquakes 

To  complement  the  above  results  from  detailed  comparison  of  four  records  with 
the  corresponding  synthetics,  a  more  statistically  oriented  analysis  was  performed  by 
comparing  the  SS/S  amplitude  ratios  from  about  80  records  with  the  amplitude  ratios 
predicted  for  the  elastic  case.  There  are  several  reasons  for  doing  this  additional 
study.  First,  long  period  S  waves  are  not  attenuated  much  in  a  high  Q  environment,  so 
measuring  the  amount  of  attenuation  is  difficult,  and  performing  many  measurements 
allows  a  more  quantitative  determination  of  statistical  bounds  on  t  *  as  opposed  to  the 
more  qualitative  results  we  were  forced  to  draw  from  the  Far-Eastern  data  above.  Also, 
it  is  difficult  to  accurately  model  the  wavetrain  of  the  large  events  needed  to  produce 
visible  multiple  S  phases  by  a  simple  double-couple  mechanism.  By  using  many  events, 
hopefully  these  variations  will  be  averaged  out  and  will  not  have  too  significant  an  effect 
on  the  attenuation  estimate.  Finally,  this  is  a  fairly  simple  way  to  estimate  the 
attenuation  and  using  this  data  set  it  is  easier  to  restrict  the  area  of  investigation  to 
the  shield  areas  proper. 

The  locations  of  the  Central  Asian  events  and  European  recording  stations  used 
are  shown  in  the  map  in  Figure  7.  The  surface  reflection  points  for  the  SS  phases  are 
also  shown  in  Figure  7;  they  are  located  in  the  central  portions  of  the  Russian  shield. 
The  source  parameters  including  the  fault  plane  solutions  are  given  in  Table  II.  When¬ 
ever  possible,  the  solutions  were  taken  from  the  literature  (references  with  Table  II). 
For  the  events  for  which  there  are  no  published  solutions,  we  have  derived  our  own 
solutions  using  P  wave  first  motions  from  the  bulletins  of  the  International  Seismic  Cen¬ 
tre.  Most  events  have  a  pair  of  complementary  focal  solutions.  In  each  case,  we  com- 


puted  the  radiation  pattern  correction  for  each  focal  solution,  and  in  nearly  all  cases, 
both  solutions  gave  similar  absolute  amplitude  ratios.  We  did  not  use  the  data  in  cases 
where  there  was  a  large  difference. 

Station-event  pairs  were  chosen  so  that  one  component  of  the  receiver  was  within 
10°  of  transverse  so  as  to  eliminate  most  P-SV  contamination  and  provide  simpler 
records.  Station-event  pairs  were  eliminated  if  the  receiver  was  close  to  a  node  in  the 
radiation  pattern  for  either  S  or  SS. 

The  double-couple  corrected  amplitude  ratios  of  SS/S  are  plotted  as  functions  of 
epicentral  distance  in  Figure  8.  The  theoretical  predictions  using  the  SNA  model  of 
Grand  and  Ilelmberger  (1984)  are  also  shown  for  comparison.  These  curves  are  the 
results  of  ray  tracing  calculations  assuming  a  lossless  medium  and  smoothed  over  3° 
intervals.  Geometrical  spreading  has  been  taken  into  account  in  the  ray  tracing  calcu¬ 
lations.  To  verify  the  validity  of  these  results  we  have  also  computed  synthetic  seismo¬ 
grams  using  full  wave  theory  (Cormier  and  Richards,  1977)  and  WKBJ  synthetics  (Chap¬ 
man,  1978;  Dey-Sarkar  and  Chapman,  1978)  (Figure  9)  for  several  distances  covering 
the  same  range  and  found  a  good  agreement  between  the  results  of  the  three  methods. 
At  distances  where  there  are  multiple  arrivals,  the  arrival  of  largest  amplitude  was 
read  off  of  the  records  and  the  largest  amplitude  arrival  was  used  from  the  ray  tracing 
calculations  or  synthetic  seismograms.  Small  perturbations  of  the  SNA  model  did  not 
change  the  overall  trend  of  these  curves,  although  in  limited  distance  ranges  the 
amplitude  ratios  varied  considerably.  Therefore,  we  have  concluded  that  our  deduc¬ 
tions  will  not  be  very  model-dependent  as  long  as  we  use  a  high  velocity,  shield  type 
upper  mantle. 

Although  there  is  a  considerable  scatter  in  the  observed  amplitude  ratios  in  Figure 
8,  the  mean  trend  of  the  data  points  is  clearly  below  the  theoretical  curve  for  t  =0 
seconds,  thus  indicating  anelastic  losses  in  the  mantle  for  SS  relative  to  S.  The  1$ 
which  fits  the  data  points  with  the  least  rms  error  is  around  2.5  to  3  seconds.  These 
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values  are  comparable,  given  the  uncertainties  associated  with  such  measurements,  to 
the  values  used  by  Burdick  et  al  (1982)  to  match  the  observed  amplitude  ratios  of  SS 
relative  to  S  in  synthetic  seismograms  for  the  same  area.  Due  to  the  averaging  over 
many  events  it  is  unlikely  that  these  results  are  biased  by  source  directivity  or  a  few 
wrong  fault  plane  solutions.  Given  this  many  events,  the  results  do  not  change  even 
when  isotropic  source  radiation  is  assumed. 

In  the  35  to  65°  distance  range,  SS  has  multiple  arrivals  due  to  interactions  with 
the  420  and  670  km  discontinuities.  For  this  study,  we  took  the  maximum  SS  ampli¬ 
tude  observed  on  each  record  in  determining  the  SS/S  amplitude  ratio.  In  most  cases, 
the  maximum  observed  amplitude  is  on  the  same  branch  of  the  SS  triplication  as  the 
maximum  SS  amplitude  on  the  synthetic  seismograms  calculated  for  an  anelastic 
earth  (£*= 0  seconds).  However,  the  branch  with  maximum  amplitude  changes  over  the 
45°  to  80°  range  as  seen  in  the  WKBJ  synthetics  in  Figure  9.  We  are  now  in  the  process 
of  recalculating  SS/S  amplitude  ratios  individually  for  the  different  branches,  and  this 
should  allow  further  refinement  of  our  t  *  models. 

We  have  also  attempted  to  estimate  t  *  from  the  decay  rates  of  amplitudes  of  suc¬ 
cessive  multiple  ScS  phases  from  Hindu-Kush  and  Far  Eastern  events  where  the  surface 
reflection  points  are  all  under  shield  areas.  The  event  parameters  are  listed  in  Table  III 
and  the  surface  reflection  points  are  shown  in  Figure  10.  Individual  amplitudes  are 
corrected  for  geometrical  spreading  and  radiation  pattern.  On  some  of  the  long  period 
records  ScS  comes  in  close  behind  S,  so  the  amplitude  ratios  were  taken  with  respect 
to  ScS2.  The  amplitude  of  the  ratio,  (amplitude  ScSN  /  amplitude  ScS2),  can  be  written 
as  A  ~  exp  [-7rf(N-2)f  *],  where  t*  is  the  t *  for  one  leg  of  the  path  (a  single  ScS).  Then, 
In  A  ~  -7Tf(N-2)f  *,  and  t  *  can  be  estimated  from  the  slope  of  a  plot  of  In  A  vs  (N-2)f,  as  in 
Figure  11.  The  t$  value  for  a  single  passage  through  the  mantle  is  ~4.2  seconds.  This 
gives  an  average  Q  of  230  over  the  path  of  a  single  ScS,  remarkably  similar  to  the  value 
estimated  by  Sipkin  and  Jordan  (1980)  for  QScS  for  continental  regions.  These  results 


deeply  enough  to  sample  such  a  layer.  It,  Is  also  of  interest  that  higher  values  of  / c 
were  found  for  paths  crossing  the  upper  mantle  under  the  tectonic  areas  of  Eurasia, 
than  under  the  shield  areas  of  Eurasia.  This  suggests  a  tg  differential  of  around  '<■’ 
seconds  for  long  period  S  waves  for  shield  versus  tectonic  paths,  corresponding  to  a  ip 
differential  of  >0.5  seconds.  Such  a  tp  differential  is  larger  than  the  similar  differential 
proposed  for  short  period  waves,  implying  t*  versus  frequency  curves  for  shield  and 
tectonic  regions  which  are  separated  over  the  0.01  to  10  Hz  band  and  divergent  at  low 
frequencies. 
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Event  Parameters 
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Figure  Captions 


Figure  1.  Event  locations,  recording  stations,  and  surface  reflection  points  of  multiple 
S  phases  for  several  deep  earthquakes  in  the  Far  East. 


Figure  2.  SH  traces  from  deep  Far  Eastern  earthquakes,  recorded  at  digital  European 
stations,  aligned  with  the  corresponding  synthetics  generated  using  full  wave  theory. 
The  synthetics  are  for  an  elastic  earth  with  the  SNA  velocity  structure  (Grand  and 
Helmberger,  1984)  with  the  appropriate  radiation  pattern  taken  into  account;  a  source 
time  function  is  not  included  in  these  synthetics.  Some  of  the  event  parameters  are 
noted  with  each  record,  and  a  fuller  description  is  in  Table  I. 


Figure  3.  Waveform  matching  to  estimate  t  *  of  SS  and  SSS  phases.  In  each  suite  of 
waveforms,  the  first  waveform  is  the  observed  phase.  The  succeeding  waveforms  are 
synthetics  which  include  a  source  transfer  function,  radiation  pattern,  and  attenuation 
per  the  value  of  t  *  listed  next  to  the  waveform.  The  t  *s  estimated  by  comparing  ampli¬ 
tudes  of  the  unattenuated  synthetics  with  the  observed  waveforms  are  a)  1  second,  b) 
1.8  seconds,  and  c)  8  seconds. 


Figure  4.  SSS/SS  spectral  amplitude  ratios  for  synthetic  waveforms  generated  for  an 
anelastic  earth.  The  spectral  ratios  are  quite  fiat,  consistent  with  the  possibility  that 
mantle  structure  makes  no  significant  contribution  to  t  *  in  the  long  period  band. 


Figure  5.  Representative  spectral  ratios  of  long  period  multiple  S  and  ScS  waves. 


Figure  6.  Examples  of  waveform  matching  for  a)  SS-SSS  and  b)  ScS-ScSg.  In  both 
cases,  the  top  trace  is  the  original  waveform,  and  the  bottom  trace  is  another  arrival 
from  the  same  record.  In  (a),  a  hilbert  transform  has  been  applied  to  SS  in  the  sefcond 
trace  to  account  for  the  phase  shift  between  SS  and  SSS.  Various  values  of  t  have 
been  applied  to  the  original  waveform  in  the  intermediary  traces,  to  find  the  value  of  t 
which  best  describes  the  broadening  between  the  two  phases. 


Figure  7.  Event  locations,  receivers,  and  surface  reflection  points  for  SS  phases  from 
central  Asian  events. 


Figure  8.  Plot  of  SS/S  amplitude  ratios  versus  epicentral  distance.  The  amplitude 
ratios  have  been  corrected  for  radiation  pattern.  The  theoretical  predictions  of  the 
SNA  model  (Grand  and  Helmberger,  1984)  for  a  number  of  t '  values  are  also  shown  for 
comparison. 


Figure  9.  Long  period  S  and  SS  WKBJ  synthetics  for  distances  between  45°  and  80°, 
generated  for  the  SNA  velocity  model  and  an  anelastic  earth.  The  ratios  of  the  max¬ 
imum  SS  amplitude  to  the  S  amplitude  are  in  good  agreement  with  the  SS/S  amplitude 
ratios  from  raytracing  in  Figure  10. 


Figure  10.  Event  locations,  recording  stations,  and  surface  reflection  points  of  multiple 
ScS  phases  with  surface  reflection  points  under  the  Eurasian  shield. 

Figure  11.  Plot  of  In  (An/A2)  versus  (N-2)f  for  ScS  phases  where  ScS,,  and  ScSN  have 
surface  reflection  points  under  a  shield  region.  The  slope  of  a  least  squat  es  fit.  through 
the  points  gives  t$  ~  4.2  seconds. 
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therefore,  in  any  study  of  frequency  dependence  of  Q,  to  restrict  the  investigation  to  a 
chosen,  Limited  area,  as  far  as  this  is  possible,  in  all  frequency  bands.  Phis  was 
attempted  in  our  study. 

Clearly,  the  specific  forms  of  this  frequency  dependence  for  various  parts  of  the 
Earth's  mantle  could  provide  valuable  information  about  the  physical  processes  occur¬ 
ring  in  the  Earth.  It  is  of  special  interest  to  find  any  relationships  between  regional 
variations  in  Q  and  regional  variations  in  velocity  and  anisotropy  such  as  those  being 
found  in  recent  tomographic  studies  of  lateral  variations  of  clastic  properties  in  the 
earth  (Clayton  and  Cromer,  1983;  Nataf  et  al,  1984;  Woodhouse  and  Dziewonski,  1984). 
Moreover,  accurate  estimation  of  yields  of  nuclear  explosions  also  requires  knowledge 
of  attenuation  properties  of  the  mantle. 

This  study  is  an  attempt  to  obtain  a  frequency  dependent  Q  model  for  the  mantle 
underlying  the  northern  shield  areas  of  Eurasia.  We  have  used  estimates  of  t  horn 
short  and  long  period  observations  over  the  0.01  to  10  Hz  band.  Details  of  the  t  meas¬ 
urements  are  in  Der  et  al  (1985)  and  Lees  et  al  (1935)  (hereafter  referred  to  as  Paper  I 
and  Paper  II,  respectively). 

Previous  work  on  the  frequency  dependence  of  Q  includes  regional  studies  across 
the  0.5  to  4  Hz  band,  especially  with  respect  to  differences  between  the  Eastern  and 
Western  United  States  (Lay  and  Helmberger,  1981;  Der  et  al,  1982a;  Der  and  Lees, 
1984).  At  longer  periods,  Sipkin  and  Jordan  (1979)  have  found  that  multiple  ScS  waves 
required  a  frequency  dependence  of  Q  in  the  western  Pacific.  Spanning  frequencies 
across  the  larger  seismic  band  of  0.01  to  10  Hz,  several  specific  frequency  and  depth 
dependent  Q  models  of  the  Earth  have  been  proposed  (Lundquist,  and  Cormier,  1980; 
Anderson  and  Given,  1982)  based  on  temperature  and  pressure  dependent  physical 
models  of  anelastic  attenuation.  The  common  feature  of  these  models  is  that  Q 
increases  with  frequency  in  most  of  the  seismic  band  and  the  depth  dependence  is 
such  that  the  upper  mantle  low  Q  zone  is  confined  to  shallower  depths  at  shorter  wave 


individual  spectral  component  is  reduced,  The  waveforms  and  spectral  ratios,  on  the 
other  hand,  are  mostly  shaped  by  t.he  apparent  t  ,  t  ,  which  for  a  limited  frequency 
band,  may  be  written  as: 

tT=t*+f(dt  7  df). 

Since  the  apparent  f.  *  is  less  affected  by  amplitude  fluctuations  due  to  random 
heterogeneities  in  the  Earth  and  is  thus  easier  to  measure,  more  apparent  f*  measure¬ 
ments  were  made  than  absolute  ones.  The  absolute  t  will  then  have  to  be  obtained  by 
using  sets  of  the  solutions  of  the  differential  equation  above  that  agree  with  all  the 
measured  apparent  t  *,  and  available  absolute  t '  as  well,  in  some  optimum  sense. 

In  using  a  variety  of  methods  we  attempted  to  reduce  the  chances  of  methodology 
related  biases  and  by  using  averages  of  measurements  we  tried  to  minimize  the  effects 
of  scatter  in  the  results  which  is  so  typical  of  all  kinds  of  geophysical  data.  Generally, 
the  observed  variance  of  the  data  dictates  the  amount  of  data  needed  to  establish 
results  with  a  desired  confidence  limit.  Unfortunately,  as  the  results  of  the  previous 
two  parts  of  this  paper  show,  the  scatter  in  the  data  is  much  worse  than  one  would  like 
it  to  be,  even  in  the  long  period  band.  A  risky  practice  prevalent  in  numerous  studies 
of  earthquake  source  mechanisms  using  band  limited  data  is  to  model  seismograms  in 
the  time  domain  by  using  complex  source  mechanisms  in  which  the  number  of  parame¬ 
ters  is  dangerously  close  to  the  degrees  of  freedom  available  from  the  band  limited 
data  themselves.  Clearly,  by  including  enough  parameters  anything  can  be  fit,  and  Q 
can  be  easily  traded  off  with  source  properties.  We  have  purposely  tried  to  avoid  this 
problem  by  canceling  the  sources  in  comparing  S  and  SS  and  multiple  ScS  amplitude 
levels  for  a  large  number  of  events,  rather  than  attempting  to  derive  detailed  source 
models  for  fewer  events.  This  way,  we  feel  that  we  made  a  cruder,  but  more  robust 
statement  about  the  amplitude  loss  in  these  phases  even  though  the  outcome  may 
appear  to  be  less  aesthetically  pleasing  to  the  reader  than  those  from  many  detailed 


waveform  studies. 


The  selection  of  the  types  of  measurements  we  have  made  was  also  dictated  by  the 
nature  of  the  data  themselves.  At  the  high  frequency  end  of  the.  spectrum,  up  to  10  Hz. 

P  waves  constitute  the  only  kind  of  teleseismie  data  available.  S  waves  usually  do  not 
eontain  detectable  energy  above  2  Hz  even  for  high  Q  paths.  At  frequencies  in  the  0.1 
to  2  Hz  range  a  variety  of  measurements  on  both  P  and  S  waves  can  be  utilized.  In  this 
band  we  may  investigate  any  deviations  from  the  assumption  tp=  4  tp  sinee  both  P  and 
S  waves  are  observable.  At  frequencies  below  0. i  Hz  it  is  not  possible  to  do  this  since 
the  amplitudes  and  speetra  of  long  period  P  waves  are  quite  insensitive  to  Q  and  any  Q 
related  changes  are  masked  by  other  effeets.  Therefore,  at  the  long  period  end  of  the 
speetrum  we  had  to  rely  exclusively  on  analyses  of  multiple  S  and  SeS  waves. 

An  important,  but  as  yet  largely  unknown  factor  in  all  studies  of  Q  is  the  effeet  of 
scattering  on  body  wave  speetra,  amplitudes,  and  waveforms.  The  various  types  of 
measurements  utilized  in  our  study  are  affected  in  various  degrees  by  scattering. 
According  to  reeent  theoretical  and  numerical  studies,  the  spectra  of  short  period 
body  waves  are  relatively  less  affected  by  scattering  than  waveforms  and  amplitudes 
(Frankel  and  Clayton,  1984;  MeLaughlin  et  al,  1984).  Thus  speetra -will  tend  to  be  more 
determined  by  anelasticity,  while  waveforms  and  amplitudes  will  be  diagnostic  of  both 
the  anelastic  and  scattering  losses.  This  seems  to  be  supported  by  the  faet  that  ana¬ 
lyses  speetral  shapes  generally  give  higher  Q  values  than  those  of  waveforms  and  wave 
amplitude.  We  may  consider,  therefore,  speetral  Q  estimates  as  upper  limits  on  Q, 
while  those  from  waveforms  (rise  times  for  instance)  will  result  in  lower  limits.  We 
believe  that  scattering  losses,  while  not  negligible,  are  smaller  than  the  anelastie  losses 
that  oecur  in  shear  deformation.  This  contention  is  supported  by  regional  and  global 
studies  where  a  4:1  ratio  between  t$  and  tp  explains  both  the  speetral  shapes  and 
amplitude  data,  while  if  scattering  were  dominant  this  would  not  be  true  (Lay  and 
Helmberger,  1981;  Der  et  al,  1982b).  A  4:1  ratio  between  tp  and  tp  also  implies  small,  if 
any,  losses  in  compression,  (Sailor  and  Dziewonski,  1978;  Anderson  and  Given,  1980). 


Figure  2  shows  the  frequency  dependent,  Q  model  for  the  Eurasian  shield,  EURS, 
that  we  developed  to  fit  the  observations  listed  above  and  shown  in  Figure  1.  Table  1 
contains  a  listing  of  the  Q  model.  Figure  3  shows  t*  and  t  predicted  by  our  Q  model 
for  direct  S  or  P  at  60°,  superposed  on  the  observations  shown  in  Figure  1.  Clearly,  the 
model  is  a  smoothed,  simple  fit  to  the  data.  Many  details  of  the  model  can  not  be  well 
resolved,  and  our  intention  was  merely  to  find  a  plausible  model  that  fits  the  available 
data.  Nevertheless,  the  resultant  model  has  some  specific  properties  that  are  required 
by  the  observational  analyses.  As  expected,  Q  generally  increases  with  frequency,  and 
it  also  increases  with  depth  for  each  frequency.  There  is  a  low  Q  layer  between  100  and 
200  km;  this  is  required  by  the  spectral  ratios  between  the  branches  of  the  short 
period  P  wave  triplications.  It  is  not  clear  whether  the  broadening  of  the  low  Q  zone  at 
100  to  200  km  is  required  by  the  long  period  data.  However,  the  upper  mantle  Q's  are 
still  high  in  an  overall  global  context.  It  is  interesting  that  the  data  seem  to  require 
this  low  Q  layer  while  the  observed  extent  of  the  triplication  branches  of  the  same  high 
frequency  P  arrivals  do  not  necessarily  require  any  kind  of  low'  velocity  layer  for  P 
waves  in  the  upper  mantle  (King  and  Calcagnile,  1976;  Paper  1).  The  EURS  model  also 
has  a  low  Q  layer  at  the  base  of  the  mantle  because  the  amplitudes  of  long  period  mul¬ 
tiple  S  arrivals  give  t$  ~  2.5  to  3j|  seconds  while  multiple  ScS  arrivals  give  ~  4.2 
seconds.  The  exact  values  of  Q  and  thickness  of  this  layer  are  pnot  well  resolved  since 
the  long  period  ScS  arrivals  tend  to  be  poor  in  continental  regions  (Sipkin  and  Jordan, 
1980). 

Figures  4  and  5  show  t§  and  tscS  ihr  the  EURS  Q  model  as  a  function  of  epiccntral 
distance  and  frequency.  If  Q  was  constant  throughout  the  earth,  t  *  would  be  directly 
proportional  to  travel  time  and  smoothly  increase  with  distance.  However,  this  work  is 
consistent  with  observations  that  t  *  is  fairly  constant  with  distance  in  the  30  to  B0Q  dis¬ 
tance  range  (Der  and  McElfresh,  1977;  Shore,  1983).  This  is  because  arrivals  at  shorter 
distances  have  traveled  more  obliquely  through  the  upper  mantle,  sampling  more  of 
the  upper  mantle  low  Q  than  arrivals  at  greater  distances  which  travel  more  steeply 
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Figure  Captions 


Figure  1.  nummary  of  t*  observations  reported  in  Der  el.  al  (1905)  and  Lees  ct  ni 
(1985).  Each  box  is  iabeled  with  the  phases  and  method  of  analysis  used.  The  sizes  of 
the  boxes  are  representative  of  the  range  of  variations  of  the  observations. 


Figure  2.  The  EURS  Q  model.  Each  line  is  a  plot  of  Q  versus  depth  for  a  different  fre¬ 
quency. 


Figure  3.  Plot  of  i  *  and  t '  for  direct  S  or  P  at  60°  as  predicted  by  the  EURS  Q  model, 
superimposed  on  the  t*  observations  of  Figure  1. 


Figure  4.  Plot  of  t  *  versus  epicentral  distance  for  direct  S  waves,  assuming  the  SNA 
shear  velocity  model  (Given  and  Helmberger,  1980)  and  the  EURS  Q  model. 


Figure  5.  Plot  of  t  *  versus  epicentral  distance  for  direct  ScS  waves,  assuming  the  SNA 
shear  velocity  model  (Given  and  Helmberger,  1980)  and  the  EURS  Q  model. 
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Abstract 


In  this  paper  we  discuss  some  aspects  of  estimating  t  *  from  short  pci  iou  body 
waves  and  present  some  limits  on  t*(f)  models  for  the  central  and  southwestern 
United  States  (CUS  and  SWUS).  We  find  that  for  short  period  data,  with  frequencies 
above  1  or  2  Hz,  while  the  average  spectral  shape  is  stable,  the  smaller  details  of  the 
spectra  are  not;  thus,  only  an  average  t* ,  and  not  a  frequency-dependent  t  ,  can  be 
derived  from  such  information.  Also,  amplitudes  are  extremely  variable  for  short 
period  data,  and  thus  a  great  deal  of  data  from  many  stations  and  azimuths  must  be 
used  when  amplitudes  are  included  in  attenuation  studies. 

The  predictions  of  three  pairs  of  models  for  t  *(/ )  in  the  central  and  southwestern 
United  States  are  compared  with  time  domain  observations  of  amplitudes  and 
waveforms  and  frequency  domain  observations  of  spectral  slopes  to  put  bounds  on  the 
attenuation  under  the  different  parts  of  the  country.  A  model  with  the  t*  values  of  the 
CUS  and  SWUS  converging  at  low  frequencies  and  differing  slightly  at  high  frequencies 
matches  the  spectral  domain  characteristics,  but  not  the  time  domain  amplitudes  and 
waveforms  of  short  period  body  waves.  A  model  with  t*  curves  converging  at  low  fre¬ 
quencies,  but  diverging  strongly  at  high  frequencies  matches  the  time  domain  observa¬ 
tions,  but  not  the  spectral  shapes.  A  model  with  nearly-parallel  t  (/ )  curves  for  the 
central  and  southwestern  United  States  satisfies  both  the  time  and  frequency  domain 
observations. 

We  conclude  that  use  of  both  time  and  frequency  domain  information  is  essential 
in  determining  t*(f)  models.  For  the  central  and  southwestern  United  States,  a  model 
with  nearly-parallel  t*(f)  curves,  where  Ai *~0.2  seconds,  satisfies  both  kinds  of  data  in 
the  0.3-2  Hz  frequency  range. 


Introduction 


In  the  last  several  years,  there  has  been  considerable  progress  in  the  study  of 
attenuation  in  the  mantle,  including  much  work  on  the  frequency  dependence  of  Q 
across  the  seismic  band.  There  is  now  a  profusion  of  frequency  dependent  t '  models  in 
the  literature  and  in  various  research  reports.  Many  of  these  have  been  derived  from  P 
wave  spectra  with  heavy  reliance  on  the  details  of  the  spectral  shapes,  especially  at  low 
frequencies.  Others  are  derived  from  broad-band  studies  of  time  domain  amplitudes 
and  waveforms  without  using  spectral  information.  A  number  of  these  studies  have  also 
shown  that  the  anelastic  attenuation  is  higher  for  paths  crossing  the  upper  mantle 
under  the  western  United  States  than  for  paths  under  the  central  and  eastern  United 
States,  and  several  pairs  of  models  for  frequency  dependent  t*  have  been  proposed  for 
such  paths. 

There  are  differences  and  inconsistencies  among  the  results  mentioned  above, 
often  due  to  interpretation  of  only  one  aspect  of  the  data  or  from  use  of  inappropriate 
methods  of  data  analysis.  Thus,  in  the  first  part  of  this  paper,  we  wish  to  point  out 
some  of  the  pitfalls  of  working  with  short  period  data  by  discussing  the  use  and  limita¬ 
tions  of  time  domain  and  spectral  domain  information  for  the  determination  of  t  *(f )  in 
the  0.3-2  Hz  range.  Short  period  time  domain  amplitude  data  is  subject  to  substantial 
variability,  so  large  amounts  of  data  must  be  averaged  over  a  range  of  azimuths  in 
order  to  draw  significant  conclusions.  The  site  dependent  fluctuations  in  the  P  wave 
spectral  shapes  are  such  that  they  preclude  the  reliable  determination  of  the  func¬ 
tional  forms  of  £*(/)  from  spectra  of  short  period  P  waves  arriving  from  a  limited 
range  of  azimuths  to  any  single  station  or  a  moderate  sized  array. 

In  the  second  part,  we  attempt  to  constrain  the  differences  in  t  *{f )  between  the 
central  United  States  (CUS)  and  the  southwestern  United  States  (SWUS)  by  examining 
spectral  and  amplitude  data  for  short  period  P  and  S  waves  in  the  two  regions.  Three 


pairs  of  models  for  the  two  regions  are  examined,  two  with  <cys( 0  and  tswvstf)  converg¬ 
ing  at  low  frequencies  and  diverging  at  different  rates  at  high  frequencies,  and  one 
where  tcusi.^)  and  ^SVVS^)  2u~c  nearly  parallel.  We  find  that  the  quasi-parallel  mouel 
best  matches  both  the  spectral  and  amplitude  data. 


Use  of  Time  and  Frequency  Domain  Information  in  D  ctcrmirung  ** 

A  variety  of  mutually  complementary  techniques  arc  available  for  determining 
t*(f )  from  such  time  and  frequency  domain  information  as  changes  in  amplitudes, 
waveforms,  and  frequency  spectra  of  body  waves.  At  a  constant  frequency,  t  can  be 
determined  from  the  change  in  amplitude  of  a  phase  (relative  to  another  phase  or  an 
assumed  source),  where 

4-g-u/r  (i) 

The  waveforms  and  spectral  ratios,  on  the  other  hand,  are  mostly  shaped  by  the 
apparent  t * ,  t',  which  for  a  limited  frequency  band,  may  be  written  as: 

F  =  t*  +  f  (dt*  /  df)  .  (2) 

F  can  be  determined  directly  from  spectral  measurements  or  spectral  ratios: 

p-  _  -1  d(ln  A)  (3) 

7 T  df  '  ■  , 

Since  the  t*(f  )  for  given  paths  in  the  Earth  probably  represent  superpositions  of 
elementary  absorption  bands  it  is  unlikely  that  they  would  correspond  in  their  func¬ 
tional  forms  to  any  simple  absorption  band  model  although  they  are  often  modeled  as 
such.  Therefore,  we  chose  to  leave  the  functional  forms  unspecified  and  to  be  solely 
determined  from  the  data. 


Short  period  body  wave  amplitudes  and  waveforms  are  characterized  by  strong 
variations  associated  with  small  changes  in  the  positions  of  the  sources  and  receivers. 


The  site  response  spectrum,  Rj,  is  responsible  for  the  strong  variations  of  P 
waveforms  and  amplitudes  between  closely  spaced  receivers  at  arrays  (Mack  19(  9, 
Frasier  &  Filson  1972;  Chang  &  von  Scggcrn  1980).  The  physical  causes  of  these  varia¬ 
tions  are  focusing  and  defoeusing  due  to  lateral  heterogeneities  in  the  Earth,  local 
conversions  to  Rayleigh  waves,  and  the  like.  Only  in  rare  instances  can  the  site  effects 
be  adequately  modeled  by  the  responses  of  plane-parallel  layering.  Since  the  nature  of 
the  site  response  function  strongly  influences  estimates  of  Q  and  thus  the  choice  of 
methodology  for  estimating  t  \  we  shall  discuss  its  properties  in  some  detail.  The  first 
important  property  of  Rj  is  that  at  most  sites  it  is  strongly  dependent  on  the  slowness 
vector  (azimuth-distance)  of  the  arrival.  The  variations  in  the  waveforms  due  to  varia¬ 
tions  in  the  slowness  vector  are  random  and  so  can  easily  be  separated  from  Q  effects. 
Another  important  property  of  the  site  response  function  is  that  In  the  spectral  domain 
its  average  slope  is  quite  stable  over  a  reasonably  wide  frequency  range  in  the  short 
period  band,  although  there  is  considerable  variation  in  the  details  of  the  individual 
spectral  shapes.  Figure  1  is  a  histogram  of  relative  t *  measurements  in  the  0.5-4  Hz 
band  between  sets  of  pairs  of  sensors  at  NORSAR  for  10  arbitrarily  chosen  seismic 
events  at  various  azimuths.  The  standard  deviation  of  this  A t  population  is  0.06 
seconds,  which  indicates  that  this  measurement,  which  reflects  the  variations  in  the 
gross  shapes  of  the  spectra,  is  quite  stable  across  NORSAR.  In  practical  terms,  it 
means  that  only  5%  of  the  measured  values  of  Af*  differ  by  more  than  0.12  seconds 
from  the  average.  This  explains  the  fact  that  apparent  t  measurements,  which 
characterize  the  broad-band  gross  spectral  shapes  in  terms  of  a  single  parameter  t* 
(per  equation  (3)),  are  usually  comparable  in  various  published  studies  (f  easier  &  Fil¬ 
son  1972;  Noponen  1975;  Der  &  McElfresh  1977).  Using  the  empirical  formula  A mb  ~ 
1.35AH*  (Der  et  a l  1979),  Af*  =  0.06  seconds  translates  into  A mb  ~  0.03.  The  actual 
variation  of  Am6  as  measured  from  amplitude  variations  is  about  0.4  magnitude  units, 
demonstrating  that  the  average  spectral  shapes  are  more  stable  than  the  signal  ampli¬ 
tudes.  Of  course,  this  experiment  does  not  rule  out  the  possibility  of  azimuthal  bias  in 


Lt  is  clear  that  within  the  variance  of  the  data  a  number  of  t *  models  can  describe  the 
spectra.  In  a  more  extensive  study,  Bache  (1984)  stacked  spectra  at  the  United  King¬ 
dom  arrays  for  suites  of  Soviet  nuclear  explosions  and  found  that  the  stacked  spectra! 
ratio  shapes  can  be  fitted  to  a  variety  of  different  tp  models  in  the  absence  of  con¬ 
straints  at  lower  frequencies.  Thus,  at  least  for  the  path  from  the  Russian  shield  to 
northern  Europe,  while  the  averaged  spectra  are  inconsistent  with  £  *  models  that  arc 
strongly  frequency  dependent  at  high  frequencies,  the  variance  in  the  spectral  shapes 
is  such  that  we  cannot  distinguish  between  a  £  *  which  is  independent  of  frequency  and 
af*  which  varies  slowly  with  frequency  in  the  1-7  Hz  band. 

The  apparent  t*  measurements,  which  are  a  reasonably  stable  estimate  of  the 
average  £ *  in  the  short  period  band,  must  thus  be  connected  in  some  way  to  other 
kinds  of  data  to  define  £*(/).  A  feasible  approach  is  to  combine  apparent  tp  measure¬ 
ments  with  short  period  S  analyses  and  any  available,  reliable  absolute  t  *  estimates  in 
the  intermediate  and  long  period  bands  (Der  et  a l  1982a).  S  waves  are  especially  suit¬ 
able  for  this  type  of  study;  observationally,  £5  ~  4£p  (Der  et  a l  1980;  Lay  &:  Helmberger 
1981;  Der  et  al  1982a),  so  for  a  given  path  length,  S  waves  are  much  more  sensitive  to 
the  effect  of  Q  than  are  P  waves.  Although  several  studies  of  Q  indicate  that  some 
losses  in  compression  may  occur  in  the  Earth  (Sailor  &  Dziewonski  1978;  Anderson  & 
Given  1982),  the  need  for  such  refinements  requires  further  observational  verification. 

It  must  be  pointed  out  that  because  of  the  focusing  effects  it  would  also  be  unrea¬ 
sonable  to  expect  a  good  correlation  between  the  apparent  £*  measured  from  spectral 
slopes,  and  the  body  wave  amplitudes  measured  at  single  sites  from  limited  source 
regions.  It  is  easy  to  imagine  scenarios  in  which  the  body  wave  arrivals  are  systemati¬ 
cally  focused  or  defocused  for  most  azimuths.  The  amplitudes  of  P  waves  are  even 
more  severely  affected  by  the  factor  Rj  than  are  the  spectral  shapes.  It  is  not  uncom¬ 
mon  to  see  systematic  variations  in  P  wave  amplitude  between  sensors  exceeding  a  fac¬ 
tor  of  4-5  across  large  arrays  for  all  arrivals  from  a  limited  source  region  (Chang  &  von 
Seggern  1980).  Again,  it  may  be  hoped  that  azimuthal  averaging  will  reduce  such 


biases.  At  any  rate,  amplitude  patterns  of  P  waves  arriving  from  a  limited  range  of 
azimuths  have  very  little,  if  anything,  to  do  with  Q  at  individual  recording  sites  (Butler 
1984).  In  addition  to  variance  due  to  focusing  and  defocusing  across  an  array,  scatter¬ 
ing  is  an  important  phenomenon  which  is  part  of  the  t  measurements,  though  it  is 
difficult  to  separate  the  scattering  and  intrinsic  components  of  attenuation  (Richards 
&  Menke  1983;  Bache  1984;  McLaughlin  &  Anderson  1984).  This  again  points  to  the  need 
to  average  amplitude  and  spectral  measurements  in  estimating  t  ,  though  there  is 
much  work  still  to  be  done  on  the  effect  of  scattering  contributions  on  t\  Neverthe¬ 
less,  it  is  possible  to  distinguish  among  various  types  of  t  (/ )  models  by  combining 
amplitude,  waveform,  and  spectral  observations  for  both  P  and  S  waves  in  the  short 
period  band  as  illustrated  below. 

Constraints  on  t*(f)  Models  for  the  Central  and  Southwestern 
United  States  using  Observations  of  Short-Period  P  and  S  Waves 

Several  pairs  of  frequency  dependent  Q  models  have  been  proposed  to  explain  the 
differences  in  observed  short  period  P  and  S  wave  amplitudes,  waveforms,  and  spectra 
between  the  central  and  southwestern  United  States.  We  shall  limit  our  following  dis¬ 
cussions  to  these  two  regions  since  they  show  the  largest  regional  differences  in  mantle 
Q  structure  in  the  United  States,  notably  larger  than  the  variations  found  between  the 
eastern  and  western  United  States.  (Solomon  &  Toksoz  1970;  Der  et  al  1975;  Butler  & 

Ruff  1980). 

These  models  fall  into  two  categories.  The  first  group  is  characterized  by  t *(/ ) 
curves  that  converge  to  the  same  value  at  low  frequencies.  Such  models  are  based  on 
the  idea  of  an  "absorption  band  shift"  where  the  low  frequency  limits  of  the  absorption 
band  are  constant  in  all  areas  and  the  decrease  of  tp  (and  t$)  with  frequency  is  deter¬ 
mined  by  a  single  parameter,  rm.  The  convergent  models  proposed  by  Lay  &  Helm- 


burger  (1981)  and  Butler  (1984)  are  nearly  identical  and  have  a  strong  frequency 
dependence  for  paths  crossing  the  mantle  under  the  eastern  United  States.  These 
models  were  developed  as  bounds  on  the  t  *  variations  in  the  United  States  to  fit  ampli¬ 
tude  data  for  frequencies  ^1  Hz.  and  are  used  as  an  example  of  a  strongly  convergent 
model.  Another  pair  of  convergent  £*(/)  models  was  proposed  by  Hadley  &  Mellman 
(1983)  to  explain  the  t  *  and  mb  differentials  between  NTS  and  RKON.  This  pair  has  a 
weak  frequency  dependence  in  order  to  fit  the  average  apparent  tp  estimates  from 
spectral  measurements  across  the  United  States. 

The  third  pair  of  f'(f)  models  do  not  converge  at  low  frequencies  and  are  quasi¬ 
parallel  throughout  the  seismic  band.  This  model  was  proposed  by  Der  t \t  al  ( 1982a, b) 
in  order  to  conform  with  some  reported  regional  Q  differentials  at  lower  frequencies 
(Solomon  &  Toksoz  1970;  Lee  &  Solomon  1975).  The  idea  of  regional  Q  differentials  over 
a  broad  frequency  range  encompassing  both  the  short  and  long  period  bands  is  further 
supported  by  regional  variations  in  the  Q  estimates  derived  from  long  period  ScS  waves 
(Sipkin  &  Jordan  1979;  Der  et  al  1984).  It  also  appears  to  be  a  more  natural  choice  if 
attenuation  is  a  thermally  activated  process,  in  which  case  temperature  increases  in 
broad  regions  of  the  mantle  would  cause  an  increase  in  attenuation  in  a  broad  fre¬ 
quency  band.  Therefore,  both  theoretical  considerations  and  the  available,  though  few, 
results  on  the  regional  variations  of  Q  in  the  long  period  band  appear  to  indicate  that 
regional  Q  differentials  will  exist  in  the  long  period  band  and  the  t*  curves  do  not  con¬ 
verge.  One  can  reach  the  same  conclusion  from  short  period  observations  alone,  as  we 
shall  show  below. 

In  Figure  4a  we  show  graphs  of  the  three  t  *(/  )  pairs  discussed  above  as  models  for 
the  frequency  dependence  of  attenuation  along  shield-to-shield  (S-S)  and  shield-to- 
tectonic  (S-T)  types  of  paths;  the  corresponding  t r  curves  are  shown  in  Figure  4b.  The 
convergent  models  with  large  and  small  t  *  differentials  at  high  frequencies  are  named 
CS  (convergent,  strong  frequency  dependence  at  high  frequencies)  and  CW  (conver¬ 
gent,  weak  frequency  dependence  at  high  frequencies),  respectively;  the  quasi-parallel 


model  is  named  QP.  Note  the  substantial  variations  in  the  0.3-2  Hz  range  of  Af  *  and  A/  * 
for  the  S-S  and  S-T  paths  of  the  different  models. 

In  this  study,  we  are  considering  the  central  areas  of  the  United  States  as  part  of 
the  "shield"-type  environment,  while  the  more  attenuating  Basin  and  Range  and  Cordil- 
lerian  areas  of  the  United  States  are  included  under  the  label  "tectonic".  The  events 
studied  are  deep  earthquakes,  so  the  source  side  of  the  raypath  is  not  subject  to  upper 
mantle  attenuation,  and  is  thus  considered  as  "shield",  since  previous  work  has  shown 
that  the  t  *  values  for  paths  originating  from  deep  earthquakes  do  not  differ  from  those 
originating  from  near-surface  events  on  shields  (Der  et  al  1982b).  We  shall  also  use 
some  results  of  regional  studies  of  P  and  S  wave  amplitude  and  spectral  variations 
between  these  two  regions  (Booth  et  at  1974;  Der  et  at  1980,  1982a;  Shore  1983;  Butler 
1984). 

Any  acceptable  pair  of  t‘(f)  models  for  shield-to-shield  and  shield-to-tectonic 
paths  must  satisfy  a  number  of  time  domain  and  spectral  constraints  derived  from 
observations.  In  this  study,  we  look  at  the  results  of  analyses  of  short  period  data  from 
WWSSN  and  LRSM  stations  in  the  United  States.  Though  the  WWSSN  and  LRSM  instru¬ 
ments  have  slightly  different  responses,  the  results  from  the  different  studies  are  simi¬ 
lar.  In  the  time  domain,  the  CUS-SWUS  contrast  is  about  a  factor  of  2  or  3  in  P-wave 
amplitudes  around  1-2  Hz  (Booth  et  at  1974;  Butler  &  Ruff  1980;  Der  et  al  1982a),  with 
a  slight  increase  of  the  average  periods  in  the  WUS  to  0.9  seconds  as  compared  to  0.7 
seconds  in  the  EUS  (Der  et  al  1982a).  S  waves  around  0.3-1  Hz  have  amplitudes  that 
average  6-10  times  larger  in  the  CUS  than  in  the  SWUS  as  seen  by  typical  short  period 
instrumentation  (Lay  &  Helmberger  1981;  Der  et  al  1982a).  For  these  S  wave  ampli¬ 
tude  ratios,  the  effects  of  sediment  amplification  are  unimportant  and  the  large 
differentials  must  thus  be  attributed  to  lateral  variations  of  mantle  Q  (Der  et  al  1982a). 
Depending  on  the  observed  S  wave  period  seen  in  the  CUS,  which  presumably  shows  the 
original  source  pulse  with  less  modification  due  to  Q,  there  is  a  varying  amount  of  pulse 
broadening  in  the  S  waves  seen  in  the  SWUS.  If  the  S  wave  has  a  high  frequency 


character  in  the  CUS  the  S  wave  period  increases  considerably  in  the  SWUS  for  thr 
same  event,  sometimes  from  1  second  to  more  than  3  seconds.  If  the  source  has  a  low 
frequency  character  no  great  change  in  the  pulse  shape  is  evident  (her  et  al  1980). 
This  is  in  agreement  with  the  results  of  numerical  experiments  on  broad-band  and 
narrow-band  signals  (Der  &  McElfresh  1980). 

In  the  spectral  domain,  estimates  of  tp  for  the  shield-to-tectonic  type  of  path  are 
around  0.45  seconds  and  for  shield-to-shield  paths  the  value  is  around  0.2  seconds  (Der 
&  McElfresh  1976,  1977;  Der  et  al  1984),  consistent  with  the  differential  in  tp  between 
the  CUS  and  the  SWUS  of  around  0.2-0.25  seconds  (Der  &c  McElfresh  1977;  Der  et  al 
1982a;  Shore  1983)  in  the  short  period  band.  The  differential  in  apparent  t$  is  about  3 
to  4  times  that  of  tp  (Der  et  al  1980,  1982a;  Lay  &  Helmberger  1981),  thus  supporting 
the  idea  that  most  differential  Q  losses  occur  in  shear  deformation.  Table  1  presents  a 
summary  or  these  time  and  frequency  domain  constraints  for  P  and  S  waves.  Figure  5 
shows  representative  S  waves  at  CUS  and  SWUS  WWSSN  stations,  with  the  appropriate 
gain  factors  noted  next  to  the  waveform.  The  difference  in  both  amplitude  and  period 
between  S  waves  from  the  same  event  arriving  in  the  CUS  and  the  SWUS  is  quite  strik¬ 
ing. 

It  is  also  interesting  to  consider  how  the  regional  variations  of  t  correlate  with 
other  geophysical  properties.  As  others  have  noted  there  seems  to  be  a  substantial 
correlation  between  the  t*  variations  discussed  above  and  the  travel  time  delays  across 
the  United  States  as  reported  by  Cleary  &  Hales  (1966),  Herrin  &  Taggert  (1968),  and 
Sengupta  &  Julian  (1976).  The  major  features  of  their  results  include  large  delays  in 
the  Basin  and  Range  and  SWUS,  early  arrivals  through  the  CUS  shield  region,  and 
slightly  late  arrivals  along  the  Atlantic  seaboard  and  in  New  England,  very  similar  to 
the  variations  in  t  *  that  were  outlined  above.  Heat  flow  and  the  available  electrical 
conductivity  data  have  also  been  found  on  a  broad  scale  to  follow  similar  patterns 
across  the  United  States  (Roy  et  al  1972;  Gough  1984).  These  correlations  of  t*  varia¬ 
tions  with  travel  time  anomalies,  heat  flow,  and  thermal  conductivity  are  consistent 


with  the  idea  that  temperature  is  an  important  parameter  in  the  physics  of  attenuation 
mechanisms. 

Now,  we  investigate  how  the  models  in  Figure  4  fit  the  variety  oT  observations  given 
in  Table  I.  Figures  6,  7,  and  B  show  simulated  waveforms  for  the  three  t  *(/ )  pairs  for  a 
variety  of  source  pulse  durations.  The  Ohnaka  pulse  (Hermann  &  Kijko  1983)  was  used 
for  the  source  function  since  this  type  of  waveform  does  not  produce  some  of  the  pecu¬ 
liar  looking  waveforms  which  can  result  from  the  corners  of  the  triangular  or  tra¬ 
pezoidal  source  pulses  that  are  commonly  used  in  calculating  such  synthetics.  The 
source  is  then  convolved  with  both  a  WWSSN  short  period  instrument  response  and  the 
appropriate  t  *(/ )  to  model  the  P  and  S  waveforms.  The  period  of  each  waveform  and 
the  amplitude  ratio  of  each  pair  of  waveforms  is  also  shown  in  the  figures. 

From  Figures  6,  7,  and  8,  it  is  clear  that  the  three  t  *(/  )  models  have  substantially 
different  effects  on  the  waveforms.  There  are  significant  differences  in  both  the  S-S/S- 
T  amplitude  ratios  and  the  amount  of  waveform  spreading  predicted  by  the  three 
models.  The  amplitude  and  period  variations  in  these  figures  were  used  to  judge  the  fit 
or  the  various  models  to  the  time  domain  constraints  outlined  above.  To  judge  the  fit  of 
the  models  to  spectral  constraints  we  use  the  apparent  t  graphs  shown  in  Figure  4b. 

The  CW  model,  with  t*  convergent  at  low  frequencies  and  small  differences  in  t*  at 
higher  frequencies  was  designed  to  fit  the  observed  t  differentials  around  1  Hz,  and 
thus  it  is  consistent  with  the  spectral  domain  constraints.  However,  the  model  fits  nei¬ 
ther  the  observed  P  or  S  wave  amplitudes  or  periods.  Because  this  model  has  high 
values  of  t '  in  the  0.5  to  2  Hz  range,  it  produces  relatively  long  period  synthetics,  even 
for  a  source  modeled  as  an  impulse.  The  shortest  period  shield-shield  P  and  S  waves 
that  this  model  predicts  are  0.75  and  2.6  seconds,  respectively,  while  much  shorter 
period  P  and  S  waves  are  observed  at  central  United  States  WWSSN  stations.  This 
model  also  predicts  S-S/S-T  amplitude  ratios  for  P  and  S  waves  of  around  1.3,  much 
less  than  observed,  especially  for  S  waves.  In  the  CW  model,  £*( f)  varies  rapidly  with 


frequency  around  1  Hz  to  give  the  desired  A t  *  differential.  Thus,  while  At  *  is  large,  A / 
is  relatively  small  as  seen  in  Figures  4a  and  b,  which  explains  the  relatively  small 

amplitude  ratios. 

The  QP  and  CS  models  are  both  consistent  with  the  observed  S-S/S-T  amplitude 
ratios  for  P  and  S  waves.  However,  a  mueh  longer  source  duration  is  needed  by  the  CS 
model  relative  to  the  QP  model  to  produce  the  appropriate  S-S  periods.  Within  the  con¬ 
siderable  scatter  of  the  data,  both  models  also  agree  with  the  difference  in  period  of  P 
and  S  arrivals  between  the  shield-shield  and  shield-teetonic  paths,  again,  for  the 
appropriate  initial  S-S  periods.  The  CS  model  tends  to  predict  a  much  larger  inerease 
in  period  between  spreading  of  pulses  whieh  arrive  in  a  teetonie  region  and  those  which 
arrive  in  a  shield  region  than  observed,  though  the  model  eannot  be  ruled  out  when  the 
variation  in  the  data  is  considered. 

The  major  differences  between  the  QP  and  CS  models  appear  when  the  spectral 
data,  or  the  relative  Ps  are  considered.  The  CS  model  prediets  AtJ  ~  0.7  seconds  or 
i,tj  ~  2.9  seconds  around  1  Hz,  about  three  times  the  observed  values.  These  very 
large  P  differentials  are  due  to  the  large  P  differentials  and  the  substantial  frequency 
dependence  around  1  Hz  in  this  model.  On  the  other  hand,  the  QP  model  prediets  A tP  ~ 
0.2  seconds  or  AtJ  ~  0.8  seconds,  in  reasonable  agreement  with  the  observations. 

While  the  QP  model  is  consistent  with  the  data,  it  can  still  be  refined.  The  t 
eurves  eould  be  somewhat  more  divergent  around  1  Hz,  slightly  increasing  the  model 
predictions  of  the  S-S/S-T  amplitude  ratios  and  of  Af  . 

This  comparison  of  the  short  period  time  and  spectral  domain  data  with  the  pred¬ 
ictions  of  different  pairs  of  P  models  emphasizes  the  importance  of  utilizing  both  the 
time  domain  and  spectral  information  in  developing  a  P  model.  The  analyses  above 
mostly  constrain  the  differentials  in  the  P(/)  for  the  two  types  of  paths,  but  a  allow  a 
considerable  variation  in  the  absolute  values  in  P(/)  at  low  frequencies.  These  values 
have  to  be  determined  independently  in  the  long  and  intermediate  period  bands. 
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Summary  and  Conclusions 


In  this  paper,  we  have  discussed  the  importance  of  utilizing  both  time  anti  spectral 
domain  information  in  developing  t  models  and  the  limitations  in  using  the  various 
kinds  of  data.  For  short  period  data,  while  the  average  spectral  shape  is  stable,  the 
small  details  of  the  spectra  are  not;  thus,  such  information  can  only  be  used  to  gen¬ 
erally  constrain  t*  and  not  to  fit  a  specific  model  to  the  data.  Also,  amplitudes  are  very 
variable  for  short  period  data  as  evidenced  from  array  recordings;  thus  to  use  short 
period  amplitude  data  in  attenuation  studies,  a  great  deal  of  data  from  many  stations 
and  azimuths  is  needed  to  have  confidence  in  the  results. 

To  put  bounds  on  the  attenuation  under  the  central  and  southwestern  United 
States,  we  compared  the  amplitude,  waveform,  and  Af  predictions  of  three  models 
with  observations  from  short  period  P  and  S  waves.  The  models  included  two  where  the 
values  of  £'  for  the  S-S  and  S-T  paths  converged  at  low  frequencies,  with  strong  and 
weak  frequency  dependence  at  the  high  frequencies,  respectively,  and  one  where  the  t* 
curves  for  the  two  paths  are  roughly  parallel.  Even  given  the  considerable  scatter  in 
the  amplitude  and  waveform  data,  the  observations  ruled  out  the  two  convergent 
models.  The  convergent  model  with  the  small  t  differential  at  high  frequencies  was 
consistent  with  observed  differentials  in  apparent  t  between  the  eastern  and  western 
United  States,  but  did  not  predict  large  enough  amplitude  or  period  variations  between 
the  EUS  and  the  WUS.  The  convergent  model  with  the  large  t  *  differential  correctly 
predicted  the  observed  amplitude  ratios  and  waveform  period  differentials  between  the 
EUS  and  WUS,  but  predicted  a  much  larger  A t*  between  the  two  regions  than  is 
observed.  The  quasi-parallel  model  provides  a  reasonable  fit  to  both  the  time  and  fre¬ 
quency  domain  observations. 

From  the  above  analyses,  we  conclude  that  models  which  converge  towards  low 
frequencies  are  not  acceptable  representations  of  attenuation  under  the  United  States, 
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Figure  1.  Relative  t  *  between  subarrays  at  NORSAR  for  ten  teleseismic  events. 


Figure  2a.  Intersensor  spectral  ratios  for  a  set  of  NTS  nuclear  explosions  at.  NORSAR. 


Figure  2b.  The  same  as  Figure  2a,  except  that  the  intersensor  spectral  ratios  are  for  a  set  ol 
Kazakh  explosions  at  NORSAR. 


Figure  3.  Array  averaged  power  spectra  at  NORSAR  of  (a)  a  Soviet  peaceful  nuclear  explosion 
on  24  November  1972  and  (b)  a  Kazakh  explosion  on  10  July  1973.  The  bounds  are  95% 
confidence  limits.  Theoretical  spectra  derived  using  a  von  Seggern-Blandford  (1972^  source 
model  and  various  t  *  models  are  superimposed  on  the  averaged  spectra.  The  other  t  models 
used  in  this  figure  are  the  quasi-parallel  (QP),  weakly  convergent  (CW),  and  strongly  conver¬ 
gent  (CS)  models,  all  for  shield-to-shield  (S-S)  paths.  These  t  *  models  are  shown  in  Figure  4a. 


Figure  4.  (a)  t*( f)  and  (b)  £*(f)  models  discussed  in  the  text  and  used  for  the  waveform  syn¬ 
thetics  in  Figures  6,  7,  and  8.  =  4 tp  has  been  assumed. 


Figure  5.  Short  period  S  waveforms  from  WWSSN  stations  in  the  United  States.  For  each 
event  noted  on  the  left,  a  waveform  is  shown  from  a  station  in  the  central  United  States 
(center)  and  a  station  in  the  southwestern  United  States  (right).  For  each  waveform,  the  sta¬ 
tion  designation,  the  short  period  instrument  orientation,  the  instrument  gain  in  thousands, 
and  the  period  of  the  S  phase  are  shown  listed  from  top  to  bottom.  The  CUS/SWUS  S  wave 
amplitude  ratio  [4s  =  ( Agz  +  A§N  +  A$B)UV\  is  noted  between  each  pair  of  waveforms  for  the 
cases  where  all  three  components  were  available  at  both  stations. 


Figure  6.  Synthetic  P  and  S  waveforms  corresponding  to  the  S-S  and  S-T  paths  for  the  QP  t 
model  shown  in  Figure  4a  and  a  short  period  WWSSN  instrument  response.  The  source  pulses 
are  shown  on  the  left;  all  of  the  time  series  are  longer  than  shown,  so  that  none  of  the  source 
or  synthetic  pulses  have  been  truncated.  The  S-S/S-T  peak-to-peak  amplitude  ratio  is  shown 
between  each  pair  of  waveforms,  and  the  period  is  written  below  each  waveform  in  italics. 


Figure  7.  Synthetic  P  and  S  waveforms  corresponding  to  the  S-S  and  S-T  paths  for  the  CS  t  * 
model  shown  in  Figure  4a  and  a.  short  period  WWSSN  instrument  response.  The  source  pulses 
are  shown  on  the  left;  all  of  the  time  series  are  longer  than  shown,  so  that  none  of  the  source 
or  synthetic  pulses  have  been  truncated.  The  S-S/S-T  peak-to-peak  amplitude  ratio  is  shown 
between  each  pair  of  waveforms,  and  the  period  is  written  below  each  waveform  in  italics. 


Figure  B.  Synthetic  P  and  S  waveforms  corresponding  to  the  S-S  and  S-T  paths  for  the  CW  t  * 
model  shown  in  Figure  4a  and  a  short  period  WWSSN  instrument  response.  The  source  pulses 
are  shown  on  the  left;  all  of  the  time  series  are  longer  than  shown,  so  that  none  of  the  source 
or  synthetic  pulses  have  been  truncated.  The  S-S/S-T  peak-to-peak  amplitude  ratio  is  shown 
between  each  pair  of  waveforms,  and  the  period  is  written  below  each  waveform  in  italics. 


